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Abstract 


A  sub-scale  nozzle/cowl  configuration  was  previously  studied  experimentally  to  determine 
the  low  speed  performance  of  a  high  speed  nozzle.  The  current  computational  study 
investigated  the  same  configuration  and  compared  the  experimental  and  computational 
data  for  the  hypersonic  nozzle/cowl  performance  at  off-design  conditions  (subsonic 
through  supersonic).  A  combination  of  Van  Leer  and  Roe  flux-splitting  algorithm, 
developed  at  WL/FIMC  was  used  to  solve  the  flow  field.  The  explicit  formulation  used 
assumed  laminar,  planar  flow  and  the  perfect  gas  equation  of  state.  A  validation  for  the 
CFD  algorithm  was  accomplished  by  comparing  the  experimental  data  according  to  the 
parametric  cases  and  the  grid  refinement.  Grid  refinement,  geometry,  and  operating 
conditions  all  influenced  the  comparison  of  the  experimental  and  computational  data. 
Generally,  as  the  Mach  number  and  NPR  (Nozzle  Pressure  Ratio)  increased,  the 
consistency  between  the  experimental  and  the  computational  results  showed  better 
agreement,  although  the  computational  investigation  was  performed  on  the  basis  of  planar, 
laminar  flow.  For  the  different  cowl  configurations,  the  results  for  baseline  and  -5  degrees 
cowl  show  a  better  agreement  than  the  long  and  +5  degree  cowl.  But  the  discrepancy 
between  the  planar  model  and  the  three-dimensional  model  caused  the  inconsistency  at 


the  secondary  separation  bubble. 


I.  INTRODUCTION 


1.1  Background 

Hypersonic  flight  is  a  vital  link  in  the  advancement  of  transportation  within  and 
beyond  the  earth's  atmosphere.  Whether  it  requires  the  relatively  short  exposure  of  a  re¬ 
entry  vehicle,  or  the  prolonged  cruise  of  a  hypersonic  transport,  surviving  the  aero- 
thermodynamic  environment  characteristic  of  this  flight  regime  constrains  the  design  of 
a  wide  class  of  aerospace  vehicles.  The  success  of  hypersonic  air-breathing  vehicles,  such 
as  the  proposed  National  AeroSpace  Plane  (NASP),  depends  largely  on  efficient  design 
of  integrated  propulsion  and  aerodynamic  systems.  Much  of  the  recent  increase  of  interest 
in  hypersonic  aerodynamics  can  be  attributed  to  the  effort  centered  around  the 
development  of  the  envisioned  NASP  (1).  The  NASP  utilizes  a  supersonic  combustion 
ramjet  (SCRAMjet)  propulsion  system.  The  location  of  this  SCRAMjet  on  the  underside 
of  the  airframe,  and  the  shaping  of  the  forebody  and  afterbody  are  of  prime  concern  to 
propulsion/airframe  integration  (2).  Figure  1-1  shows  an  integrated  airframe-propulsion 
concept  for  a  hypersonic  airplane,  wherein  the  entire  undersurface  of  the  vehicle  is  part 
of  the  SCRAMjet  engine. 

Integrating  airframe  and  propulsion  systems  in  transonic  and  supersonic  flight  has 
been  an  area  of  active  research  over  the  past  decade.  Because  these  flowfields  are 
complex  and  defy  solution  by  simple  analytical  methods,  the  preponderance  of  reports 
represent  joint  experimental/computational  studies.  As  discussed  in  Reference  3,  a  vehicle 
such  as  NASP,  which  would  be  required  to  operate  efficiently  at  high  speed,  would  be 


heavily  influenced  by  the  aero-thermo  design  constraints  of  hypersonic  flight.  However, 
with  the  additional  requirement  of  operating  from  a  runway,  the  subsonic  and  transonic 
regimes  of  its  flight  envelope  also  become  important  considerations  (1). 

The  hypersonic  propulsion  system  operates  at  extremely  high  pressure  ratios 
(SCRAMjet  total  pressure  divided  by  ambient  pressure,  PTJPJ  at  hypersonic  speeds.  The 
nozzle/afterbody  geometry  is  very  important  in  efficiently  extracting  thrust  from  the 
propulsion  system.  For  hypersonic  speeds,  the  exhaust  flow  is  generally  underexpanded, 
and  exerts  high  pressure  on  the  afterbody,  which  translates  into  additional  thrust.  At 
subsonic  and  low  supersonic  Mach  numbers,  the  nozzle/afterbody  is  operating  at  an 
extremely  overexpanded  condition,  and  the  pressures  in  the  afterbody  are  below  ambient 
pressure,  causing  drag  and  reducing  effective  thrust  (4). 

A  major  element  of  the  integrated  design  is  the  interaction  between  the  multiple 
jet  plume  flow,  the  external  flow  and  the  afterbody.  Unfortunately,  major  portions  of  the 
hypersonic  flight  envelope  cannot  be  simulated  by  existing  ground-test  facilities. 
Therefore,  accurate  and  validated  computational  fluid  dynamics  tCFD)  codes  are  needed 
to  provide  simulations  of  realistic  flight  conditions.  The  experimental  validation  data 
should,  of  course,  model  as  many  of  the  realistic  flight  conditions  as  possible.  The 
validated  codes  can  then  be  used  to  provide  the  most  reliable  predictions  of  the 
increments  in  performance  or  design  parameters  associated  with  the  differences  between 
the  available  test  conditions  and  the  flight  environment  ( 1 ).  This  research  addresses  the 
nozzle/afterbody  integration  of  the  hypersonic  propulsion  system  at  low  speed,  off-design 
conditions  between  Mach  numbers  of  0.6  to  3.0. 
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1.2  Purpose 


Propulsive  nozzles  for  hypersonic  vehicles  represent  an  extremely  demanding  test 
for  CFD  codes.  The  llowfield  of  a  propulsive  nozzle  is  a  complicated  structure  with 
expansion  waves,  shock  wa\es.  contact  surfaces,  and  the  interaction  among  all  three  and 
surface  boundaries  such  as  the  nozzle  wall  or  cowl.  Several  CFD  codes  have  been 
developed  to  resolve  the  propulsive  nozzle  llowfield  of  a  hypersonic  vehicle. 

One  such  CFD  algorithm  used  to  obtain  a  computational  flow  solution  is  an 
explicit  flux-vector  splitting  (FVS)  and  flux-difference  splitting  (FDS)  finite  volume 
scheme.  The  FVS  algorithm  is  quite  robust  near  sonic  points  with  continuously 
differentiable  flux  contributions  and  the  FDS  (Roe-averaged  scheme)  is  particularly 
desirable  due  to  its  shock  capturing  capability  and  stability  at  high  Mach  numbers  (5). 

The  purpose  of  this  research  is  to  compare  results  between  experimental  and  CFD 
data  for  hypersonic  nozzle/cowl  performance  at  off-design  conditions  (subsonic  through 
supersonic).  The  development  of  the  NASP  will  be  based  on  the  effective  combination 
of  CFD  and  experimental  testing  to  arrive  at  the  aircraft's  final  configuration  efficiently 
and  with  a  large  degree  of  confidence. 

1.3  Scope 

This  research  effort  compares  the  computational  data  obtained  from  the  flux 
splitting  code  with  experimental  data  for  a  hypersonic  propulsive  nozzle  flowfield. 
Cochran  (4)  conducted  experiments  investigating  the  low  speed,  off-design  performance 
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of  a  generic  two-dimensional  hypersonic  nozzle/afterbody  configuration  in  the  Wright 
Laboratory  (WL)  Two-Foot  Trisonic  Gasdynamic  Facility. 


In  the  numerical  method,  the  flow  in  the  two-dimensional  hypersonic  nozzle  is 
assumed  to  be  unsteady  and  laminar  at  zero  degrees  angle  of  attach  and  may  have 
complex  shock  wave  interactions.  It  also  involves  strong  shock  -  laminar  boundary-layer 
interactions  which  may  result  in  separated  regions.  The  major  interests  in  the  CFD 
simulation  data  of  the  experimental  results  will  be  the  comparison  of  the  length  and 
thickness  of  the  separation  bubble  on  the  top  surface  of  the  nozzle.  The  fluid  in  the  nozzle 
was  treated  as  a  calorically  and  thermally  perfect  gas. 

The  present  investigation  of  the  nozzle/cowl  CFD  code  was  conducted  for  three 
parametric  cases:  Mach  number  effects,  nozzle  cowl  configuration  effects,  and  nozzle 
pressure  ratio  (NPR)  effects.  The  parametric  test  matrices  for  numerical  implementation 
are  presented  in  Table  1-1.  The  Mach  number  effects  have  a  Mach  number  range  of 
subsonic  and  supersonic.  The  subsonic  cases  included  Mach  numbers  of  0.6  and  0.8, 
while  the  supersonic  cases  were  Mach  numbers  of  1.9  and  3.0.  To  determine  the  effects 
of  the  nozzle  lower  cowl,  four  different  cowl  configurations  were  tested:  the  baseline  (BL, 
no  extended  cowl)  (Figure  1-3),  the  long  (L,  0.5  inch  straight  extension)  (Figure  1-4),  the 
positive  5  degrees  (+5,  0.5  inch  extension  angled  5  degrees  toward  the  nozzle  centerline) 
(Figure  1-5),  and  the  negative  5  degree  (-5,  0.5  inch  extension  angled  5  degrees  away 
from  the  nozzle  centerline i  (Figure  1-6)  Photographs  of  the  cowls  installed  on  the  model 
are  shown  in  Figure  1-2.  The  NPR  effects  for  the  baseline  cowl  configuration  were 
conducted  for  Mach  numbers  0.6.  0.8,  1.9,  and  3.0.  respectively.  The  NPR  is  the  ratio  of 
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nozzle  exhaust  flow  total  pressure  (PTo)  and  the  tunnel  static  pressure. 

The  initial  external  and  internal  (nozzle)  flowfield  values  were  determined  by 
applying  the  uniform  freestream  conditions  on  the  nozzle/afterbody  ramp  surface  and  the 
internal  nozzle  flow  conditions  at  a  five  probe  total  pressure  rake  (Figure  1-7).  The  rake 
location  is  3.25  inches  upstream  of  the  nozzle  throat  to  measure  the  internal  nozzle  flow 
total  pressure.  The  nozzle  flow  total  pressure  equals  the  average  of  these  five  probes. 
Nozzle/afterbody  ramp  surface  pressure  were  measured  at  desired  NPR.  The  significant 
area  of  interest  in  comparison  of  numerical  and  experimental  model  was  the  consistency. 
The  algorithm,  implemented  by  Gaitonde  (7),  is  also  quite  robust  relative  to  grid  skewness 
and  is  accurate  to  second  order  in  both  space  and  time.  Numerical  pressure  from  the  CFD 
code  will  be  compared  to  the  experimental  data.  Additionally,  Mach  numbers  will  be 
obtained  by  the  numerical  solutions.  The  calculations  and  post  processing  were  performed 
on  the  Silicon  Graphics  Workstations  at  the  AFIT  Computer  Lab. 

1.4  Approach 

The  performance  of  two  flux-split  algorithms  in  unsteady  two-dimensional  perfect 
gas  laminar  flows  are  investigated  for  this  research.  The  flux- vector  split  method  of  Van 
Leer  and  the  flux-difference  split  method  of  Roe  are  applied  due  to  their  demonstrated 
accuracy  in  the  prediction  of  similar  high  speed  flows. 

The  flowfield  was  initially  calculated  with  the  Van  Leer  method,  which  exhibits 
shock  structures  spanning  no  more  than  two  interior  cells,  whereas  Roe's  scheme  may 
violate  the  entropy  condition  near  sonic  or  stagnation  points.  Numerical  convergence  is 
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based  on  the  computed  pressure  and  heat  transfer  on  nozzle  surface.  When  the  solutions 
in  the  investigation  properly  converged  to  a  steady  state,  the  control  parameter  for  flux 
splitting  method  is  converted  from  Van  Leer's  scheme  to  Roe’s  flux-difference  splitting, 
which  is  based  upon  accurate  prediction  of  wave  interactions  between  interfaces  through 
a  linearized  equation. 

A  sequence  of  grids  is  utilized  to  perform  a  grid  resolution  study  and  the  results 
are  compared  with  experimental  values.  The  formulation  is  finite-volume  with  second 
order  accuracy  obtained  with  MUSCL  (Monotonic  Upstream  Schemes  for  Conservation 
Laws)  approach  in  conjunction  with  a  limiter  to  prevent  oscillations.  The  advantages  of 
the  MUSCL  approach  usually  accrue  for  two  reasons  (5).  The  fluxes  are  split  according 
to  the  local  ceil  interface  Mach  number.  Also,  the  split  fluxes  are  generally  less 
differentiable  than  conserved  variables  when  transitioning  through  sonic  and  stagnation 
points.  Viscous  terms  are  centrally  differenced:  the  solution  is  marched  in  time  to  steady 
state  with  a  two  stage  Runge-Kutta  scheme.  An  implicit  Gauss-Seidel  line  relaxation 
scheme  is  also  incorporated  in  the  code  (7),  although  this  scheme  was  not  used  in  the 
present  investigation. 

The  centerline  pressures  on  the  afterbody  ramp  were  used  for  data  comparisons 
between  numerical  method  and  experimental,  because  of  the  higher  density  of  available 
data  in  this  region  (4).  The  main  parameter  for  comparison  was  the  pressure  coefficient 
of  the  nozzle/afterbody  ramp,  Cp ,  in  dimensionless  form,  which  was  calculated  from  the 
locally  measured  static  pressure,  the  upstream  wind  tunnel  test  section,  and  the  dynamic 
pressure,  respectively.  Numerical  results  were  obtained  using  the  finite  volume  scheme. 
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where  the  flow  properties  indicate  values  at  each  cell  center. 
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Figure  1-2  Hypersonic  Nozzle/ Afterbody  Model  Installed  in  the  Trisonic  Gasdynamics 


Facility 
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Table  l-l  Parametric  Cases  for  Numerical  Methods 


NPR  0.9 _ MACH  0.6  MACH  0.8  MACH  1.9  MACH  3.0 

3.0 _ 06 _ 08 _ L9 _ 10 

5.0 _ 06 _ 08 _ L9 _ 10 

7.0 _ x _ x _ L9 _ 10 

12.0  x  x  1.9  3.0 


MACH  0.6 

BASELINE 

LONG 

+5  DEG 

-5  DEG 

0.8 

BASELINE 

LONG 

+5  DEG 

-5  DEG 

1.9 

BASELINE 

LONG 

+5  DEG 

-5  DEG 

3.0 

BASELINE 

LONG 

+5  DEG 

X 

09 

_L9 

3.0 


0.9 


10 

10 

5.1 


5.0 

7.0 

12.2 


6.1 

12.1 

20.3 


II.  GOVERNING  EQUATIONS  AND  NUMERICAL  METHODS 
2.1  Governing  Equations 

The  compressible  form  of  the  Navier-Stokes  equations  adequately  describes 
aerodynamic  flow  at  standard  temperature  and  pressure.  By  neglecting  body  forces  and 
heat  sources,  the  unsteady,  viscous  compressible  form  of  the  Navier-Stokes  equations 
may  be  written  (10): 


i£  .££“u0 

dt  dXj 

“i  +  dp  lift  +  _  dry  =  Q 

dt  dXj  dxt  dXj 


(Continuity  equation ) 


(Momentum  equation) 


de  d(Et  +  p)  aCujty  -  q.) 

—  +  — -  -  — — J-  =  o  (Energy  equation ) 

dt  dXj  ctj 

P  =f(e,  P)  =/(—  -  p)  (Equation  of  state) 

P  2 

where  indicial  notation  has  been  used  (repeated  indices  are  summed)  and 
h,  =  component  of  velocity  in  the  x,  direction 
Xh}  =  viscous  stress 


(2-1) 

(2-2) 

(2-3) 

(2-4) 


duk  du{ 

T*  =  ^  +  % 


dx, 


where  5,  f  =  the  Kronecker  delta. 

ji,  X  =  first  and  second  coefficients  of  viscosity,  respectively. 


(2-5) 
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qj  =  -k(dT/dXj) 


where  k  =  coefficient  of  thermal  conductivity 

The  two-dimensional  Navier-Stokes  equations  in  Cartesian  coordinates  can  be 
written  in  vector  form  as 


W+§i+d<i=  o  (2-6) 

dt  dx  dy 


where  the  conserved  variables  are  0  =  { p,  pu,  pv,  £,}T  and  the  flux  vectors  are  separated 


into  inviscid  (F,  G )  and  viscous  (F„  ,  Gv)  terms. 


The  inviscid  flux  vectors  are 


pu 

pv 

pu2+p 

pwv 

puv 

G  = 

pv2+p 

(  E+p  )u 

(  E+p  )v 

and  the  viscous  fluxes  are 


0 

0 

Ev  = 

V 

<?v  = 

w*f 

,dT 

l^V^j 

(2-7) 


(2-8) 


The  perfect  gas  equation-of-state  closes  the  set  of  equations 

p  =  (  y-l)pe  e  =  Et  /p  -  |(«2+y2)  (2-9> 
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The  molecular  dynamic  viscosity  ji  is  given  by  Sutherland's  law  (6): 


H  =  C, 


3 

_iL 

r+c. 


(2-10) 


where  Ct  and  C:  are  constants.  For  air,  these  values  are  C,  =  2.27E-08  lbm  /(ft-sec.Rl/2) 
and  C,=  198.9  R. 


The  conservation  law  form  of  the  Navier-Stokes  equations  given  in  Eqi  itions  (2-6) 
to  (2-8)  is  more  flexible  than  that  in  Equations  (2-1)  to  (2-3)  for  purposes  of  coding  the 
Navier-Stokes  equations  into  a  CFD  algorithm:  especially  for  the  finite  volume 
applications. 

Transforming  the  Navier-Stokes  equations  from  the  physical  space  to  the 
computational  coordinate  system  is  necessary  (11).  This  transformation  drastically 
simplifies  the  applications  of  the  boundary  conditions  and  may  include  various  options 
on  grid  point  clustering  and  orthogonality,  both  being  extremely  important  for  the  solution 
of  the  Navier-Stokes  equations,  so  that  the  equations  are  consistent  with  the  computational 
domain.  Equation  (2-6)  can  be  transformed  for  use  in  a  general  coordinate  system  (£,  q) 
as  follows  (6): 


dt  0§  dr\ 


(2-12) 
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where 


u  =  rxu 


(2-13) 


F  =  J-'ftJf-FJ  * 

(2-14) 

0  -  ♦  rt,(G-GJ] 

(2-15) 

and  J  is  the  Jacobian  of  the  transformation: 


=  Vn  ~  Ve  (2'l6) 

As  a  result,  the  general  transformed  form  of  Navier-Stokes  equation  for  two-dimensions 
is  (12): 


It)  *  i£w  - ^ *  5>(G  -  G’>1 

*  -  F->  *  ">(G  -  G>)l| =  0 


(2-17) 


2.2  Discretization 

The  code  developed  at  WL/FIMC  and  applied  in  the  present  research  uses  several 
methods  in  conjunction  with  each  other  to  discretize  the  Navier-Stokes  equations.  The 
viscous  terms  in  the  tlux  vectors  are  simply  handled  by  central  differences.  The  inviscid 
terms  are  handled  by  a  combination  of  flux  vector  splitting  and  flux  difference  splitting 
as  described  below. 

Consider  the  one-dimensional,  inviscid  model  equation: 


(2-18) 


dU  +  BFiU)  = 
dt  dx 


A  first-order  explicit  discretitization  of  Equation  (2-18)  with  the  forward  Euler  method 
produces 


U 


/l+l 


vr 


F\  -  F\ 

i+-  i-- 

2  2 


=  0 


(2-19) 


Af  Lx 

where  FlW/C  is  derived  from  some  combination  of  state  of  the  flow  at  adjacent  grid  points 
[eg.  Ft^in  =  V2(F,  +  F1+i  )].  Rearranging  Equation  (2-19)  for  the  explicit  discretization 


yields 


pn  —  pn 

(t>U%  _  4  4  (2-20) 

A  t  Ax 


where 


fit/"  =  {/"+1  -  Un 


(2-21) 


The  two-dimensional  system  in  Equation  (2-6)  is  discretized  in  a  similar  manner. 
In  order  to  compute  the  right  hand  side,  or  residual,  of  the  two-dimensional  version  of 
Equation  (2-20),  the  fluxes  in  each  direction  are  successively  balanced.  Thus 


F;+1(t/")  -  F.  i (U") 


LHS, 


LHS2  =  LHS,  + 


Ax 

Ghx(U”)  '  "> 


(2-22) 


Ay 


where  LHS,  refers  to  the  left  hand  side  of  the  two-dimensional  version  of  Equation  ( 2-20). 
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The  left  hand  side  of  Equation  (2-20)  is  evaluated  by  a  two-stage  Runge-Kutta 
scheme  (Heun's  method),  while  the  right  hand  side  of  Equation  (2-20)  is  evaluated  by  a 
combination  of  Van  Leer's  flux-vector  splitting  and  Roe's  flux-dilference  splitting 

methods. 

The  two-stage  time  stepping  scheme  is  iormulated  as  (13): 

5  Un  =  '  ~  (F\  -F“)  (2-23) 

A* 

in  which 

=  j7" 

C/C— i,D  =  un  -  a, A t  (F* f  -  F^ f  ) 

*  2  2  (2-24) 

u*+ 1-2)  =  u*  -  aaA t  (F^  -  Fci}  ) 

i+2  *~2 

yn+l  _  fj(.n+l2) 

where  the  second  term  in  the  superscript  denotes  each  stage  and  its  coethcients  are  a,- 
1.0,  ou  =  0.0  . 

For  robustness,  especially  in  supersonic  expansion  regions,  the  right  hand  side  ol 
Equation  (2-20)  is  initially  calculated  using  Van  Leer's  flux-vector  splitting  method. 
Equation  (2-20),  where  the  inviscid  flux  is  split  into  positive  and  negative  components  lor 
appropriate  upwind  differencing: 

F  i  =  F*(Ul)  *  F'iU*)  (2-25) 

(+- 
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When  the  solution  properly  converged  to  the  steady  state,  the  flux  splitting  scheme  is 
convened  to  Roe's  scheme  ( 14): 


f(ul )  +  f(ur)  -  - A{u*  -  uL) 

2 


(2-26) 


The  terms  UR  and  UL,  vectors  of  conserved  variables  which  describe  the  state  of  the  flow 
at  each  interface,  are  obtained  from  the  MUSCL  approach  described  in  Section  2.3.2.  The 
"A"  indicates  a  Roe  averaged  term,  also  described  in  Section  2.3.2. 


2.3  Numerical  Flux  Split  Scheme 

2.3.1  Flux-Vector  Splitting  Scheme  (Van  Leer) 

Flux  vector  splitting  is  a  widely  used  technique  not  only  to  make  computational 
schemes  more  robust,  but  also  to  improve  the  computational  efficiency  of  finite-difference 
schemes  with  varying  degrees  of  success.  The  basic  aim  is  to  split  the  flux  vectors  such 
that  an  upwind  finite-difference  scheme  may  be  used  at  all  points  within  the  flow.  This 
is  performed  by  the  simple  method  of  separating  the  flux  vectors  into  upwind  and 
downwind  parts. 

A  description  of  upwinding  may  be  made  by  reference  to  the  following  inviscid 
equation  (5): 

+  §£.  +  $9.  =  0  (2-27) 

dt  dx  dy 

where  F  and  G  are  the  inviscid  flux  vectors.  Denoting  the  flux  Jacobians  by  A  and  B 
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respectively,  where  for  the  first-order  homogeneous  property,  F  =  AU  and  G  =  BU, 


i.e. 


F = AU ; A 


dF 

dU 


G  =  BU  ;  B  = 


BG 

dU 


A  finite-difference  form  of  Equation  (2-27)  can  be  written  as: 


At 


DFi  DGf 
Ax  Ay 


(2-28) 


(2-29) 


where  5  and  D  are  the  time  and  space  difference  operators,  respectively: 


-4jl 


(2-30) 


In  the  Steger  and  Warming  scheme,  the  first-order  homogeneous,  hyperbolic 


property  of  the  inviscid  tlux  is  utilized  to  split  the  flux  Jacobian  A, 


A  =  QlAQ  (2-31) 

where  A  is  a  diagonal  matrix  consisting  of  the  eigenvalues  of  A,  A  =  diag{«,  u+a,  u-a). 
The  eigenvalues  can  be  decomposed  into  non-negative  and  non-positive  components 


X,  =  Xf  ♦  A,-  <2-32> 

where 

X*  =  Xi  *  N  (2-33) 

2 

similarly,  the  eigenvalue  matrix  A  can  be  decomposed  into 
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A  =  A+  +  A- 


(2-34) 


where  A*  is  made  of  the  non-negative  contributions  \*  and  A'  is  constructed  of  the  non¬ 
positive  contributions  This  splitting  of  the  eigenvalue  matrix  allows  the  flux  vector 
F  to  be  rewritten  as 


FJ+1  =  F\Ul)  +  F~{UR)  =  Al+Ul  +  A**U*  (2-35) 

2 

For  second  order  accuracy,  the  vector  of  conserved  variables,  U  R  and  U  L,  are 
obtained  at  the  cell  interfaces  by  extrapolation  to  the  cell  surface  with  the  MUSCL 
approach  of  Van  Leer.  In  order  to  avoid  oscillations  in  shock  regions,  a  limiter  function 
that  employs  local  conditions  is  used  as  follows  (8): 


UL  =  U,  +  -A.  i  ; 
‘  2  l~2 


A.  l  =  g(A  i  ,  A  0  (2-36) 


u*  ■  u,.i  -  \hi  i 

4  2 


3  =«(A  3,  A  0  (2-37) 

2  *  2  J  2 


where  g  is  the  limiter  and  A,w/>  =  Ut+i  -  V, .  For  example,  Van  Leer's  limiter  is 


8  = 


A.  i  A.  i  + 

i+_  i—i 
2  2 


2  2 


A.  i  +  A.  i  +  e 
2  2 


(2-38) 


where  £  is  a  small  number  to  prevent  division  by  zero.  With  the  addition  of  the  minmod 
limiter,  the  algorithm  reverts  to  first  order  accuracy  at  shocks  in  order  to  preserve 
monotonicity  within  the  solution. 

Note  that  for  locally  supersonic  and  sonic  flow  on  both  sides  of  the  interface  (5) 
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(2-39) 


FL*  =FL  .  fr  -  =  q  if  Mi  >  !  . 

Fi+=0  ;  F*'  =  Fr  if  Ml,M*<-  1 

where  M  represents  the  local  Mach  number.  This  method  is  quite  dissipative  and  also 
exhibits  a  flux  discontinuity  where  the  eigenvalues  change  sign. 

One  cause  of  inaccuracy  in  the  Steger  and  Warming  splitting  is  that  the  split-flux 
derivatives  are  discontinuous  when  any  of  the  eigenvalues  A.,  (A)  changes  sign.  This  gives 
rise  to  instabilities  in  the  numerical  solution,  unless  some  extra  smoothing  is  built  into  the 
flux  formula.  To  solve  this  problem.  Van  Leer  developed  a  continuously  differentiable 
splitting  technique. 

Van  Leer's  approach  was  to  split  the  fluxes  so  that  the  forward  and  backward  flux 
contributions  transitioned  smoothly  at  eigenvalue  sign  changes,  i.e.  near  sonic  and 
stagnation  points.  Just  as  for  the  Steger  and  Warming  splitting,  it  was  required  that  the 
Jacobian  matrices  ( dFVdU  )  have  non-negative  eigenvalues  and  (dF'/dU  )  have  non¬ 
positive  eigenvalues  so  that  upwind  differencing  could  be  used  for  spacial  derivatives.  In 
addition,  it  was  required  that  both  Jacobians  have  one  zero  eigenvalue  for  I  M  I  <  1, 
which  leads  to  steady  transonic  shock  structures  with  only  two  interior  zones. 

The  formulas  for  F~  in  one  dimension  are  locally  given  in  terms  of  the  Mach 
number:  for  supersonic  and  sonic  flow,  I  M  I  >  1,  we  again  have; 


♦ 

11 

• 

Fr  " 

=  0 
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and  for  subsonic  flow,  1  M  1  <  1 

(2-40) 
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/f[(Y-D«  ±  2a]  I  y 

/Tv 

yf[{(Y-l)«±2a}2  /  {2(y2-1)}+v2/2] 

where 

/*  =  ±pa  [— (M±l)]2  (2-42) 

2 

2.3.2  Roe's  Flux-Difference  Splitting  Scheme 

The  Roe  scheme  for  a  one  dimensional  system  was  presented  in  Equation  (2-26), 

F",  -  [F*1  +  K  -  *  |S|  -  Op]  (2-43) 

2 

Once  again,  higher  order  accuracy  and  limiting  are  achieved  with  the  MUSCL 
approach  in  conjunction  with  a  limiter.  Yee  (15)  presents  the  MUSCL  scheme  using  a 
simple,  one  dimensional  hyperbolic  equation.  Starting  from  Equation  (2-43),  the  MUSCL 
scheme  replaces  the  UlW  and  t/,  terms  with  U  Rl+l/2  and  UL,+IP_  ,  respectively.  UR  and 
UL  are  defined  for  second-  and  third-order  spatial  differencing  by  using  various  slope 
limiters  to  eliminate  unwanted  oscillations.  A  widely  used  limiter  is  the  minmod  limiter, 
which  modifies  the  upwind-biased  interpolation  as  follows: 

V*l  =  UM  -  7  fd-<l)A,.l  +  d+i)v|  (2-44) 

*+2  4  2  2. 


(2-41) 
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UL 1  =  u,  -  — 

4  1  4 


(M)Vj.  +  d+Tj)A.+1 
2  2 


The  spatial  order  of  accuracy  is  determined  by  the  value  of  rj: 
fj  =  -1  :  fully  upwind  scheme 
r\  =  0  :  Fromm  scheme 
i*|  =  1/3  :  third-order  upwind-biased  scheme 
f]  =  1  :  three-point  central-difference  scheme 


(2-45) 


A  i  =  minmod  (A  i  ,  to  AO 


(2-46) 


A  i  =  minmod  (A  \  ,  coA  3) 
2  J*2  J+2 


(2-47) 


where 


and 


minmod  fotoy)  =  sgn(x )  .  max  {  0  ,  mini;  |x| ,  coy  JgnOO] }  (2-48) 


1  5  to  <; 


tj  *  1 


(2-49) 


Note  that  sgn(x)  means  the  sign  of  the  variable  x. 

The  symbol  "A"  in  the  above  equations  refers  to  Roe  averaging,  which  is  given  by: 
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£  _  y/Pi  +  \[Pr  xr 

y S  + 


x  =  p,  «,  v,  and  A 


(2-50) 


The  subscripts  in  Equation  (2-50).  L  and  R.  refer  to  the  components  of  the  like 
subscripted  vector  U  in  Equations  (2-44)  and  (2-45). 
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HI.  ANALYSIS  OF  NOZZLE/COWL  FLOW  FIELDS 
3.1  Nozzle/cowl  Model 

The  experimental  model  under  consideration  was  analyzed  by  Cochran  in  1988, 
who  had  as  his  goal  the  measurement  of  the  nozzle/afterbody  pressure  drag  force.  This 
force  was  calculated  from  the  integrated  external  ramp  pressure  coefficients  (Cp)  which 
resulted  from  the  impingement  of  an  overexpanded  exhaust  plume  of  the  flow  downstream 
of  the  nozzle  throat. 

The  generic  hypersonic  nozzle/afterbody  was  attached  to  a  square  forebody  (cross 
section)  with  a  15  degree  half-angle  conical  section.  Figure  1-2  is  a  photograph  of  the 
model  installed  in  the  Trisonic  Gasdynamics  Facility.  The  conical  nose  section  had  a 
sharp  nosepiece  for  supersonic  testing  and  a  round  nosepiece  for  subsonic  testing.  The 
model  lengths  were  23.33  and  22.16  inches  for  the  supersonic  and  subsonic 
configurations,  respectively. 

The  two-dimensional  hypersonic  nozzle/afterbody  was  designed  to  be 
representative  of  the  external  geometry  of  a  hypersonic  vehicle.  Figure  3-1  is  a  side-view 
drawing  of  the  nozzle/afterbody  with  the  pertinent  dimensions.  The  nozzle/afterbody  had 
a  straight  20  degree  expansions  ramp  and  a  nozzle  expansion  area  ratio  of  approximately 
6.2,  The  internal  nozzle  was  convergent  and  to  facilitate  fabrication  resulted  in  a  50 
degree  straight  convergent  ramp.  The  nozzle/afterbody  had  a  sharp  comer  at  the  nozzle 
throat  to  make  the  performance  independent  of  throat  radius.  To  determine  the  effects  of 
the  nozzle  lower  cowl,  four  different  cowl  configurations  were  tested:  the  baseline,  the 
long,  the  positive  5  degrees,  and  the  negative  5  degrees.  The  cowls  have  expansion  area 
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Figure  3-1  Hypersonic  Nozzle/ Afterbody  Dimensions  and  Coordinate  System 
(Physical  Domain) 
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ratios  of  1.0,  1.4,  1.3,  and  1.5,  respectively.  The  cowls  were  0.125  inches  thick  with 
squared-off,  blunt  trailing  edges. 

The  tests  were  conducted  in  the  WL  two-foot  trisonic  gasdynamics  tunnel 
operating  at  subsonic,  transonic,  and  supersonic  speeds.  The  maximum  Reynolds  number 
was  2.5  million  per  foot  for  the  subsonic  speed  and  5.0  million  per  foot  for  the  supersonic 
speed,  and  the  maximum  dynamic  pressure  was  350  pounds  per  square  foot  (psf)  for  the 
subsonic  speed  and  1,000  psf  for  the  supersonic  speed.  The  stagnation  temperature  was 
maintained  at  approximately  100  degrees  F.  Pressure  taps  were  placed  in  a  grid  pattern 
downstream  of  the  nozzle  assembly  in  the  internal  convergent  ramp  and  the  afterbody 


Figure  3-2  Hypersonic  Nozzle/Atterbody  Tap  Locations 
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expansion  ramp.  In  all,  sixteen  pressure  taps  were  located  along  the  centerline  of  the 
external  ramp.  Tap  locations  are  shown  in  Figure  3-2  for  the  nozzle/afterbody  left  and 
right  sides,  respectively. 

Ramp  pressure  distributions  were  obtained  over  a  range  of  NPR  for  each  Mach 


Figure  3-3  Computational  Domain  of  Nozzle/ Afterbody 


number,  cowl  configuration,  and  Reynolds  number.  The  range  of  NPR  varies  for  the 
subsonic  and  supersonic  Mach  numbers  because  of  the  constant  pressure  air  supply  for 
the  nozzle  exhaust  flow  and  the  variation  of  tunnel  static  pressure  with  Mach  numbers. 
For  each  test  configuration,  the  tunnel  flow  conditions  (Mach  and  Reynolds 
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numbers)  were  established,  and  data  for  zero  nozzle  exhaust  flow  were  obtained.  The  high 
pressure  air  was  then  established  through  the  model  and  nozzle  exhaust  flow  set  to  the 
desired  NPR.  Data  were  obtained  at  each  NPR.  For  each  test  point  the  tunnel  and  nozzle 
flow  parameters  were  recorded  during  a  scanivalve  reading  and  used  tor  the  data 
reduction. 

The  main  performance  parameter  for  the  experimental  teat  was  the 
nozzle/afterbodv  pressure  drag  coefficient  which  was  calculated  from  the  integrated 
external  ramp  pressure  coelticients.  C?,  by  assigning  an  appropriate  area  to  each  pressure 
tap  location  and  resolving  the  area  into  the  projected  axial  component.  The  centerline 
pressure  were  used  because  of  the  higher  density  ot  pressure  taps  which  characterize  the 
nozzle  exhaust  flow  and  the  afterbody  ramp  pressure  more  accurately.  The  centerline 
pressures  were  also  used  for  the  test  parameters  (Mach  number,  cowl,  and  NPR  data 
comparison). 

The  maximum  errors  measured  in  the  experimental  tests  were  0.207c  for  freestream 
total  pressuro,  0.127c  for  freestream  static  pressure.  0.57%  for  local  pressure  at  probe,  and 
0.43%  for  total  pressure  of  nozzle  exhaust,  respectively. 

3.2  Numerical  Problem  Formulation 

3.2.1  Geometry  Definition  and  Grid  Generation 

The  parametric  variations  performed  in  the  experiment  made  it  possible  to  access 
the  effects  of  nozzle  afterbody,  Mach  number,  nozzle  pressure  ratio,  and  nozzle/cowl 
geometry  on  the  forces  and  moments  induced  by  the  impinging  exhaust  plumes. 
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Modelling  the  same  effects  computationally  presented  an  opportunity  to  test  the  numerical 
method  on  a  very  practical  application.  The  geometry  for  the  nozzle/afterbody  lent  itself 
naturally  to  a  Cartesian  grid  topology,  and  the  two-dimensional  grid  was  generated  by 
GR1DGEN,  a  interactive  graphics  based  method  developed  by  the  Fort  Worth  Division 
of  General  Dynamics  Corporation  for  the  Aeromechanics  Division  of  the  Flight  Dynamics 
Laboratory  ,  Wright  Research  and  Development  Center  (16),  now  WL/FIM. 

GRIDGEN  is  made  up  of  four  modules  -  GRIDBLOCK.  GRIDGEN2D, 
GRIDGEN3D,  and  GRIDVUE3D.  These  modules  enable  users  to  create  three-dimensional 
multiple  block  structures  grids.  However,  for  this  problem,  a  quasi  two-dimensional 
hypersonic  nozzle/afterbody,  only  a  single  two-dimensional  surface  grid  is  required, 
GRIDGEN2D  was  all  that  was  required  to  develop  this  surface  grid. 

GRIDGEN2D  is  an  interactive  program  written  primarily  in  FORTRAN  and  C  for 
use  on  Silicon  Graphics,  Inc.  (SGI)  3000  and  4D  series  IRIS  Workstations  and  IBM  RISC 
6000  workstations.  In  GRIDGEN2D,  the  distribution  of  points  on  each  of  the  four  edge 
boundaries  of  a  surface  is  performed,  resulting  in  the  generation  of  surface  grids  via 
algebraic  and  elliptic  grid  generation  techniques.  The  block  dimension  (^max ,  q^  ,  Cmu) 
is  entered  to  run  in  a  single  block.  Since  a  single  surface  grid  is  desired,  =  1  was 
entered  to  eliminate  the  third  dimension.  For  the  primary  mesh,  dimensions  of  101  x  71 
x  1  were  set  in  the  %  ,  q,  and  C  directions,  respectively. 

Figure  3-3  shows  the  gnd  information  pertinent  to  the  inviscid  (Euler)  and  viscous 
(Navier  Stokes)  perfect-gas  CFD  simulations.  As  mentioned  in  (17),  dividing  the  spatial 
gnd  into  two  or  more  zones  tor  subfaces)  offers  several  advantages  for  solving  complex 


flow  problems.  For  example,  it  permits  the  domain  to  be  segmented  into  logical  subunits. 
Appropriate  placement  of  zones  also  eliminates  some  grid  singularities,  such  as 
discontinuous  body  cross  sections.  Finer  grids  can  be  used  in  multiple  zone  calculations 
than  in  single  zone  calculations  if  the  computer's  main  memory  limits  the  problem  size, 
because  only  one  zone  resides  in  main  memory  at  a  time.  It  is  assumed  that  the 
nozzle/cowl  geometry  is  composed  of  four  main  zones:  the  internal  nozzle  flow,  the 
external  flow,  the  solid  body,  and  the  plume  (jet)  boundary  zone.  The  major  section  of 
concern  is  the  internal  nozzle  flow  zone,  which  is  divided  into  two  subfaces:  surface  (4) 
and  surface  (5).  The  boundary  between  surfaces  (4)  and  (5)  (i.e.  the  throat)  was  made 
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straight  in  an  attempt  to  capture  the  sonic  condition.  At  the  solid  body  surface,  the 
appropriate  boundary  conditions  were  applied. 

The  generation  of  grid  points  on  a  surface  begins  with  the  distribution  of  points 
on  the  four  bounding  edges  of  the  surface.  When  the  coordinates  (,v,y)  have  been  set  on 
all  four  edges,  interpolation  methods  are  employed  to  distribute  points  on  the  surface 
interior. 

There  are  three  main  steps  involved  in  the  edge  generation  process.  The  first  step 
is  to  define  the  shape  of  the  nozzle/afterbody  geometry  on  which  the  edge  grid  points  are 
to  be  distributed.  Only  straight  line  segments  were  needed  to  define  the  shape  of  the 
nozzle/afterbody.  In  the  second  stage,  each  edge  shape  is  divided  into  a  number  of 
subedges,  separated  by  selected  datapoints.  called  breakpoints,  according  to  the 
experimental  configuration.  By  breaking  the  edge  into  subedges,  greater  control  of  grid 
point  placement  on  the  edges  was  gained.  In  the  third  stage,  grid  points  are  distributed 
along  each  subedge. 

Next,  algebraic  methods  are  used  to  define  grid  points  in  the  interior  of  each 
surface.  Among  types  of  algebraic  solvers,  Transfinite  Interpolation  (TFI)  with  arclength- 
based  interpolants  was  used.  Essentially,  this  method  uses  interpolants  based  on  the 
relative  arclength  between  edge  points,  which  means  that  points  close  together  on  the 
edge  remain  close  together  on  the  interior  of  the  surface.  This  method  worked  very  well 
for  the  majority  of  cases. 

Further  refinement  of  the  surface  interior  points  is  made  by  accessing  the  elliptic 
PDE  solvers.  On  surfaces  (4)  and  (5),  GRIDGEN2D’s  elliptic  solver  was  used  to  reduce 
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grid  skewness  near  the  throat  area.  Thomas-Middlecoff  eontrol  functions,  with  a  relaxation 
factor  of  0.3.  were  used  for  the  smoothing,  because  they  provided  the  best  results.  The 
small  relaxation  factor  was  chosen  so  that  the  changes  to  the  grid  would  be  slow  enough 
that  when  the  grid  was  acceptably  smooth  the  process  could  be  stopped.  Approximately 
200  iterations  were  required  to  produce  an  acceptable  grid.  The  grid  points  were  clustered 
not  only  near  body  and  surface  to  provide  accurate  resolution  of  the  boundary  layer  for 
viscous  calculations,  but  also  near  the  region  where  the  shock  is  expected  to  form. 

Figure  3-3  shows  perspective  views  of  the  grid  topology,  which  were  segmented 
into  five  subfaces  for  a  baseline  cowl  and  six  subfaces  for  long,  +5  degree,  and  -5  degree 
cowls.  Subfaces  (4)  and  (5)  had  51x36,  where  51  points  in  the  streamwise  (<;)  direction, 
36  points  in  the  normal  (r|)  direction  to  the  stream.  For  a  baseline  cowl,  subface  (2)  had 
51x10.  But  subface  (2)  had  11x10,  subface  (3)  had  41x10  for  long,  +5  degree,  and  -5 
degree  cowls,  respectively.  Subface  (1)  had  101x27.  The  grid  totaled  7,171  points.  Figure 
3-5  to  Figure  3-8,  depicting  the  grids  on  the  different  nozzle/cowl  configurations,  are  the 
output  of  GRIDGEN2D,  with  ghost  points  (points  within  the  body  surface)  added  as 
required  for  the  computational  solver.  By  using  a  preprocessor,  convgrd.f,  the  grid  points 
are  converted  into  the  cell  centered  points.  Further,  this  also  inverts  the  nozzle/afterbody 
configuration  for  subsequent  convenient  examination  of  developing  phenomena  with  the 
Flow  Analysis  Software  Toolkit  ( FAST)  package  under  development  at  NASA  Ames 
Research  Center. 

The  ghost  points  mentioned  above  are  required  because  cell-centered  finite  volume 
codes,  such  as  the  one  used  here,  solve  for  ilux  as  across  cell  faces,  which  requires 
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knowing  the  flow  variables  at  cell  centers,  rather  than  at  the  grid  points.  Ghost  points  are 
added  to  all  boundaries,  so  that  cell  centers  can  be  calculated.  Figure  3-4  illustrates  the 
addition  of  ghost  points  to  a  simple  grid,  and  the  resulting  network  of  cell  centers. 


3.2.2  Boundary  and  Initial  Conditions 

For  each  flow,  the  boundaries  may  be  categorized  into  the  following:  inflow 
boundaries,  surface  boundaries,  and  outflow  boundary. 

At  inflow  boundaries,  the  flow  vector  {p,  pu.  pv,  pe}  is  determined  according  to 
the  given  external  and  internal  freestream  values.  The  freestream  conditions  were  imposed 
on  the  cmia  and  the  faces.  Fluxes  on  the  freestream  boundaries  were  calculated  using 
the  values  for  Mach  number,  temperature,  and  pressure  assigned  in  the  initial  conditions. 
On  the  face,  the  inflow  conditions  are  specified  by 


-  U- 

where  V  is  the  vector  given  in  Equation  (3-1). 

Similarly,  the  freestream  condition  is  given  on  the  qinax  face  by 


V, 


UJmax 


=  Um 


(3-1) 


(3-2) 


Alternatively,  the  zero  gradient  condition  V,  jmax+l  =  U^jmax  was  used. 

The  boundary  conditions  at  the  nozzle/cowl  surface  are  more  complicated.  The 
viscous  terms  are  calculated  directly  from  the  conditions  specified  at  the  surface  for 
temperature  and  zero  velocity.  Since  the  computational  grid  is  cell  centered,  the  conditions 
at  the  surface  are  imposed  by  setting  the  conditions  at  the  ghost  points  just  within  the 
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surface. 
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boundary,  thus 


U, 


lmax*l,J 


=  u, 


(3-7) 


Initial  conditions  are  normally  prescribed  for  each  set  of  calculations  by  assuming 
the  freestream  conditions  exist  at  all  the  grid  points  except  at  the  boundaries  where  proper 
boundary  conditions  are  applied.  Initial  properties  -  freestream  pressure,  Mach  number, 
and  static  temperature  -  for  the  external  and  internal  (nozzle)  flow  were  taken  from 
experimental  data  at  the  rake  location  which  was  3.25  inches  upstream  of  the  nozzle 
throat.  All  properties  except  the  nozzle  static  pressure  were  assumed  to  be  constant, 
whereas  the  static  nozzle  pressure  was  assumed  by  taking  the  averaging  value.  Also,  the 
nozzle  static  top  and  bottom  pressure,  which  was  within  about  5 %  error  compared  to  the 
average  value,  was  calculated  to  verify  the  significance  of  the  initial  conditions. 


3.3  Computer  Code  Description 

Wright  Laboratory's  two-dimensional,  flux-splitting, Navier-Stokescode,  developed 
at  WL/FIMC,  was  basically  used  to  solve  the  viscous,  hypersonic,  nozzle/cowl  problem 
based  on  the  laminar  flow  assumption.  This  code  was  written  in  FORTRAN  and  was  run 
on  a  Silicon  Graphics  IRIS  Workstations. 

Data  file,  cnldat,  was  required  for  execution  of  the  code.  The  cnldat  contains 
information,  such  as  the  version  of  the  code  to  run  (Navier-Stokes  or  Euler),  whether  to 
use  an  implicit  or  an  explicit  formulation,  initial  conditions,  boundary  conditions,  and  so 
forth.  The  cnldat  file  has  essentially  four  parts:  solution  integration  parameters,  flowfield 
conditions,  boundary  conditions,  and  control  parameters.  A  sample  cnldat  file  can  be 


found  in  Appendix  A. 

Table  3-1  lists  the  parameters  that  control  the  implementation  of  the  solution, 
particularly  such  things  as  which  flux  split  scheme  is  used,  i.e.  Van  Leer,  Steger  and 
Warming,  or  Roe.  Other  parameters  include  the  ending  iteration  number,  whether  the  flow 
model  is  Euler  or  Navier-Stokes,  what  the  CFL  stability  criteria  is,  and  how  it  is  handled. 

Table  3-2  details  the  flow  field  conditions,  particularly  the  external  freestream 
conditions  and  the  internal  nozzle  freestream  conditions.  Table  3-3  lists  the  boundary 
condition  parameters,  and  the  node  range  in  which  they  apply.  Table  3-4  lists  the  parame¬ 
ters  which  control  the  iterations  between  printing  convergence  data,  the  grid  formatting 
from  ASCII  to  binary,  and  the  output  file  formatting  to  three-dimensions  to  monitor  the 
developing  situation  of  the  simulation. 


Table  3-1,  Solution  Integration  Parameters 


ISWVL 

Finite  Difference  Scheme:  3= Van  Leer,  4=Roe 

NEND 

Number  of  iterations 

INS 

Flow  model:  0=Euler,  l=Navier-Stokes 

IL,  JL 

Grid  Dimensions 

CFLMAX 

Maximum  CFL  number  allowed 

CFL 

Starting  CFL  number 

CFLEXP 

Number  of  iterations  before  CFL  doubles 

ICFL 

Number  of  iterations  between  CFL  increases 

IMPLT 

Implicit  vs.  Explicit:  0=Explicit,  l=Implicit 
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Table  3-2,  Flowfield  Conditions 


IADBWL 

Flag  to  indicate  adiabatic  wall  BC 

TWALL 

Wall  temperature 

ALPHA 

Freestream  Angle  of  attack 

RM 

Freestream  Mach  number 

REL 

Freestream  Reynolds  number  per  foot 

TINE 

Freestream  Temperature 

ALNLPHA 

Nozzle  Angle  of  Attack 

RINM 

Nozzle  Mach  number 

RINEL 

Nozzle  Reynolds  number  per  foot 

TIN  INF 

Nozzle  Temperature 

Table  3-3,  Boundary  Conditions 


INS 

Start  and  end  points  for  a  c  Solid  Body  BC 

JNS 

Start  and  end  points  for  an  rt  Solid  Body  BC 

Table  3-4,  Control  Parameters 


IREAD 

()=Deadstart.  l=Restart 

IGR1D 

Grid  format 

1P3DOP 

Plot3D  output  file  format 

MODPR 

Iterations  between  printing  convergence  data 

3.4  Solution  Procedure 

The  calculations  to  be  described  were  all  obtained  with  the  explicit  option  of  the 
code.  The  numerical  method  used  in  this  research  is  second-order  accurate  in  space  and 
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second-order  accurate  for  time.  Convergence  is  enhanced  through  the  use  of  local  time¬ 
stepping  methods. 

This  algorithm  provide  two  important  advantages.  First,  the  use  of  upwind  flux 
approximations  for  the  inviscid  components  (including  pressure)  of  Navier-Stokes 
equations  provide  a  method  for  numerically  capturing  shocks  without  the  need  for 
additional  artificial  damping  terms  which  require  adjustment  by  the  user.  In  addition,  an 
explicit  integration  scheme  provides  an  extremely  efficient  numerical  method  on  vector 
or  parallel  machines  for  solving  systems  of  equations,  because  the  dependent  variables  can 
be  explicitly  updated  using  concurrent  machine  operations. 

The  conservative  variables  at  interior  grid  points  must  be  initialized  to  start  the 
iterative  solution  algorithm  1 17  j.  The  initial  condition  for  flow  variables  was  implemented 
by  using  their  freestream  properties.  Since  this  condition  is  far  from  the  converged  result, 
a  large  perturbation  occurs  at  the  first  time  step.  Stability  bounds  constrained  the  time  step 
to  a  small  value  initially.  It  was  then  gradually  increased  as  the  solution  progressed  to  the 
steady  state.  The  code  used  in  this  investigation  allows  specification  of  the  Courant- 
Friednchs  and  Lewy  (CFL)  number  ( 19),  defined  as 


CFL  = 


Af 


(3-8) 


where  is  the  largest  eigenvalue  of  the  flux  Jacobians  matrices.  The  time  step  is  then 
locally  calculated  from  the  prescribed  CFL  number  and  the  maximum  eigenvalue.  The 
CFL  number  was  initially  set  to  a  value  0.01,  and  increased  to  0.9  for  the  maximum  value 
using  the  explicit  version  of  the  code. 
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Two  criteria  determined  convergence  of  the  solution:  the  integrated  root  mean 
square  {RMS)  pressure  and  heat  transfer  values  over  the  nozzle  surface  are  monitored. 
Both  must  remain  constant  at  convergence  with  P  and  Q  denoting  pressure  and  heat 
transfer,  respectively. 

«?W  =  jj- 

<P)*“  '  li 

where  IL  is  the  number  of  points  in  £,  direction. 

Convergence  is  assumed  when  {P)RMS,  (Q)Ri1s  do  not  fluctuate  more  than  0.05%  over  a 
large  number  of  iterations  (typically  1.000  at  CFL  of  0.9.  local  time  stepping). 

The  solution  of  the  nozzle/cowl  flowfield  started  from  freestream  conditions  is 
initiated  by  Van  Leer's  flux-splitting  with  continuously  differentiable  flux  contributions 
at  sonic  or  stagnation  point  up  to  1,000  iterations.  Tnall  was  assumed  to  be  400.6  R.  After 
that,  the  solver  was  switched  to  Roe's  scheme  which  predicted  wave  interactions 
accurately,  for  up  to  2,000  iterations.  Solutions  were  run  on  the  Silicon  Graphics 
Workstations. 

The  pressure  coefficients  were  calculated  from  locally  calculated  static  pressure 
for  comparison,  using  the  following  relationship 


where  p„  and  are  the  tunnel  test  section  static  and  dynamic  pressure,  respectively. 
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3-5  Baseline  Cowl  Grid 


Figure  3-6  Long  Cowl  Grid 
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Figure  3-7  +5  Degree  Cowl  Grid 
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IV.  COMPUTATIONAL  RESULTS 

Navier-Stokes  solutions  based  on  an  explicit  flux-splitting  scheme  were  obtained 
for  the  quasi  two-dimensional  hypersonic  nozzle/afterbody  at  freestream  Mach  numbers 
of  0.6.  0.9.  1.9,  and  3.0  at  an  angle  of  attack  of  zero  degrees.  The  external  flow,  the  jet 
exhaust  plume,  and  the  internal  nozzle  flow  were  included  in  the  calculations.  Nozzle 
Pressure  Ratios  (NPRs)  of  0.8  to  20.0  were  investigated.  In  addition,  four  different  cowl 
configurations  were  tested.  In  the  calculations,  laminar  flow  was  assumed  over  the  entire 
external  surface  and  internal  nozzle. 

As  a  necessary  step  to  the  nozzle/afterbody  investigation,  the  convergence  of  the 
numerical  solution  is  based  on  the  integrated  RMS  pressure  and  heat  transfer  over  the 
nozzle  surface.  Additionally,  grid  independence  of  numerical  results  was  tested  with  a 
series  of  grids.  A  grid  refinement  was  performed  at  M„=1.9,  NPR=7.15  and  M^,=3.0, 
NPR=16.0  on  the  long  cowl.  The  fine  grid  had  dimensions  of  201  points  in  the 
streamwise  (^)  direction,  and  151  points  in  the  normal  (r|)  direction  to  the  stream,  for  an 
overall  total  of  30,351  points.  Next,  the  medium  grid  had  dimensions  of  151  (4)  and  101 
(r|),  for  a  total  of  15,251  points.  Finally,  the  coarse  grid  had  dimensions  of  101  (<;)  and 
71  (r|)  for  a  total  of  7.171  points.  The  fine  and  medium  grids  based  on  the  cell  center 
for  the  finite  volume  scheme  are  shown  in  Figures  4- 1  and  4-2. 

The  solution  of  the  nozzle/cowl  flowfield  was  initially  performed  with  Van  Leer's 
flux  splitting,  for  the  first  1000  iterations  (Figure  4-3).  Fortunately,  solutions  in  the 
present  investigation  converged  to  a  steady  state,  using  local  time  stepping,  within  this 
1000  iterations.  Finally,  the  solution  was  fully  converged  with  Roe’s  scheme.  Solutions 
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were  performed  on  the  Silicon  Graphics  IRIS  Workstations;  the  code  expends  about  0.23 
microseconds  of  Central  Process  Unit  (CPU)  time  per  iteration  per  grid  point  for  the 
baseline  cowl  and  0.27  microseconds  for  the  long  cowl  (A/j=1.9,  NPR=7.0).  This 
calculation  required  1.216  hours  of  CPU  for  Van  Leer  and  1.233  hours  for  Roe,  totalling 
2.449  hours. 

As  expected,  the  hypersonic  nozzle  exhaust  llow  was  dominated  by  an  initial 
overexpansion  of  the  flow  downstream  of  the  nozzle  throat  for  all  test  conditions. 
Overexpansion  occurs  when  the  atmospheric  or  ambient  pressure  is  greater  than  the 
internal  exhaust  pressure,  causing  flow  to  disturb  within  the  nozzle.  The  sudden  increase 
in  the  exhaust  flow  area  of  the  divergent  afterbody  expansion  ramp  and  the  effect  of  the 
sharp  comer  at  the  throat  caused  the  flow  to  overexpand  to  the  low  pressure  and  moved 
the  separation  bubble  near  the  throat.  The  recompression  characteristic,  after  this 
expansion,  depended  on  the  parametric  cases;,  Mach  number,  cowl  configuration,  and 
NPR. 

Having  gained  some  appreciation  for  the  dynamics  of  overexpanded  jet 
impingement,  for  the  supersonic  cases,  it  is  now  of  interest  to  investigate  grid  resolution 
effects,  influence  of  the  initial  conditions,  and  parametric  variations  of  nozzle  operating 
conditions.  But  for  the  subsonic  cases,  there  was  inconsistency  between  the  experimental 
and  the  computational  results  due  to  discrepancy  in  the  models,  computational  planar 
model  and  experimental  three-dimensional  model.  Therefore,  the  explanation  for  the 
computational  results  is  focused  on  the  supersonic  cases  at  this  section.  As  parameters 
change,  the  plume/afterbody  interaction  will  be  altered,  impacting  overall  performance. 
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4.1  Grid  Resolution  Effects 

To  study  the  effects  of  grid  refinement  on  numerical  results  and  the  convergence 
of  the  numerical  solution,  a  grid  refinement  was  performed  on  a  series  of  grids,  101  x  71, 
151  x  101,  and  201  x  151..  at  M„=1.9.  NPR=7.0  and  M„=3.0,  NPR=16.0  on  the  long 
cowl.  The  long  cowl  was  chosen  because  the  numerical  solution  is  well  consistent  for  the 
baseline  cowl  with  the  assumption  of  laminar  flow.  Any  finite-dilference  method  or  finite 
volume  carries  with  it's  truncation  errors  arising  from  approximation  to  the  differentiations 
1 17).  Obtaining  numerical  solutions  on  a  senes  of  grids  tests  the  accuracy  and  uniqueness 
of  the  results  on  the  basis  of  laminar  boundary  layer. 

Comparison  of  pressure  coefficients  among  solutions  from  the  three  grids  and 
experimental  data  are  shown  in  Figure  4-4  (M„=l.9,  NPR=1.0).  The  coarse  grid  (101  x 
71)  and  medium  grid  (151  x  101)  results  follow  similar  trends  as  the  experimental  data, 
but  fail  to  place  the  reexpansion  before  the  peak  of  the  compression.  However,  the  fine 
gnd  (solid  line  in  Figure  4-4)  showed  that  the  exhaust  flow  initially  expanded  and 
compressed  downstream  of  the  sharp  throat.  Also,  it  showed  that  the  flow  re-expanded 
slightly  and  recompressed  until  the  flow  reached  the  apex  of  the  separation  bubble,  whose 
trends  were  not  revealed  in  the  coarse  and  medium  grids.  For  the  second  case  (M„=3.0, 
iVjP/?=l6.0),  the  gnd  refinement  results  showed  the  same  trends  as  the  previous  numerical 
results  (Figure  4-5).  With  the  encouraging  results  above,  if  the  finer  grid  can  be 
generated  m  the  numencal  implementation,  the  results  should  provide  better  agreement 
with  the  expenmental  data  and  more  accurately  analyze  the  development  of  the  flowfield 
phenomena.  But  even  though  the  finer  grid  was  generated  in  this  research,  there  were  still 
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difference  between  the  results  at  the  separation  bubble  due  to  the  difference  in  the 
assumption,  which  is  the  laminar  flow  in  the  numerical  performance  and  the  turbulent 
flow  on  the  experiment.  The  reference  (18)  showed  favorable  results  between  the  laminar 
and  turbulent  tlow  under  the  same  conditions,  i  e.,  two-dimensional  computational  model. 
Therefore,  it  was  determined  that  the  laminar  How  assumption  used  in  the  CFD  results 
was  the  dominant  factor  for  the  initial  separation  bubble  regime.  The  remainder  of  the 
research  was  then  performed  using  the  medium  and  fine  grid. 

The  convergence  for  the  coarse  grid,  101  x  71,  was  already  mentioned  above. 
The  medium  grid  converged  to  the  steady  state  near  2,000  iterations  and  the  fine  grid 
required  up  to  5,000  iterations.  The  code  requires  appioximately  0.14  microseconds  of 
CPU  time  per  iteration  per  grid  point  for  the  medium  grid  and  0.04  microseconds  for  the 
fine  grid  on  the  long  cowl  t.W„=3.0.  NPR=  12.0).  This  calculation  required  5.439  hours 
of  CPU  time  for  medium  grid  and  23.261  hours  for  the  fine  grid  on  the  Silicon  Graphics 
IRIS  Workstations. 

4.2  Influence  of  Initial  Conditions 

To  examine  the  influence  of  the  initial  conditions  in  the  numerical  implementation, 
different  values  for  nozzle  static  pressure  -  the  average,  top,  and  bottom  values  -  were 
used  considering  the  two-dimensional  effects  only.  With  regards  to  the  tolerances  of  errors 
from  the  experimental  data,  pressure  was  within  the  approximately  57c  and  temperature 
was  constant. 

For  this  purpose  the  NPR  and  Reynolds  number  varied  along  with  the  initial 


conditions  similar  to  tolerances  for  the  experimental  data.  As  with  the  grid  resolution 
effects,  the  long  cowl  was  chosen,  since  the  numerical  solution  was  well  consistent  to  the 
experimental  results  for  the  baseline  cowl.  Pressure  coefficient  comparing  among  solutions 
in  the  different  initial  conditions  shown  in  Figure  4-6  (M„=1.9,  NPR=1.0).  Obviously,  it 
is  very  hard  to  distinguish  the  differences  between  the  different  initial  conditions  on  the 
numerical  solutions.  These  results  indicate  that  varying  the  initial  conditions  within 
experimental  error  doesn't  significantly  effect  the  results  in  the  numerical  solutions.  A 
further  study  of  very  different  initial  conditions  may  have  produced  better  agreement 
between  the  CFD  and  experimental  data,  but  that  was  not  the  purpose  of  this 
investigation. 

4.3  Mach  Number  Effects 

The  Mach  number  effects  on  the  pressure  coefficient  distributions  are  presented 
in  Figures  4-7  to  Figure  4- 1 1  (Cp  versus  model  station  plots)  at  specific  NPRs  for  the 
baseline  cowl  configuration.  Mach  number  and  pressure  distribution  contours  over  all  the 
nozzle/afterbody  area  are  presented  in  Appendices  C  and  D. 

Data  for  the  supersonic  conditions  (Mach  number  1.9  and  3.0).  Figure  4-7  (NPR 
approximately  equal  to  one,  i.e.,  flow  off),  shows  that  the  external  flow  overexpands  to 
the  negative  Cp  then  compresses  at  the  same  location  by  the  separation  bubble 
encountered  downstream  of  the  sharp  throat.  The  data  at  Mach  number  1 .9  shows  a  higher 
rate  of  overexpansion  and  compression  than  for  Mach  3.0  condition.  For  these  supersonic 
conditions,  the  pressure  distributions  appear  nearly  flat  after  compression. 
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For  the  NPR  of  3.0  and  5.0,  Figures  4-8  and  4-9  show,  for  the  supersonic 
conditions,  that  the  exhaust  flow  continues  to  overexpand  to  a  negative  Cp  with  even 
higher  rate  than  NPR  1.0  conditions.  The  Mach  number  1.9  flow  expands  to  a  lower  Cp 
than  Mach  number  3.0,  but  continues  to  have  a  higher  rate  of  compression.  The  Mach  1.9 
Cp distribution  shows  some  recompression  characteristics  near  .t=5.0  caused  by  the  second 
separation  bubble,  but  the  exhaust  flow  does  not  recover  fully  to  the  ambient  pressure. 
The  NPR  5.0  flow  has  a  higher  rate  of  compression  than  the  NPR  3.0.  The  Mach  3.0  Cp 
distributions  are  nearly  flat. 

For  Figures  4-10  and  4-1  i,  nozzle  pressure  ratios  of  7.0  and  12.0,  respectively,  the 
pressure  distribution  characteristics  for  these  supersonic  Mach  numbers  show  similar 
tendencies  to  the  supersonic  pressure  distributions  discussed  above.  For  the  Mach  number 
1.9  pressure  distributions,  both  figures  show  an  initial  compression  of  the  exhaust  How 
downstream  of  the  throat.  The  exhaust  flow  again  reexpands  to  a  slightly  lower  pressure, 
but  not  to  the  ambient  pressure.  The  pressure  distribution  for  the  Mach  number  3.0 
condition  (Figure  4-10)  shows  the  same  compression  downstream  of  the  throat,  the  slight 
reexpansion  and  the  recompression.  The  NPR  12.0  flow  shows  a  higher  rate  of 
compression  than  the  NPR  7.0  for  both  case. 

Generally,  as  the  NPR  increases,  the  exhaust  flow  downstream  of  the  throat  has 
the  tendency  to  compress  at  a  higher  rate,  that  is,  to  form  the  thicker  separation  bubble 
due  to  the  sharp  corner  at  the  throat  for  a  given  Mach  number  conditions.  Especially  the 
exhaust  flow  reaches  to  above  ambient  pressure  for  the  NPR=12.0.  In  addition,  the  flow 
is  affected  by  the  secondary  separation  bubble  causing  the  flow  to  recompress.  Finally 
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the  nozzle  exhaust  flow  for  the  supersonic  Mach  numbers  appears  to  be  uniform. 
Thereafter,  the  flow  does  not  continue  to  expand  to  lower  pressure. 

Pressure  and  Mach  number  contour  plots  for  Mach  numbers  1.9  and  3.0,  are 
presented  in  Appendices  C  and  D.  for  the  baseline  cowl  at  the  assumed  NPR  schedule. 
In  addition  to  the  nozzle  exhaust  flow  re-expansion  and  recompression  region,  the  plume 
shock  (external  exhaust  shock),  the  plume  flow  boundary,  and  the  shocks  in  the  jet  plume 
flow  can  be  seen. 

4.4  Nozzle  Cowl  Effects 

The  nozzle  cowl  configuration  effects  for  the  supersonic  case  are  presented  in 
Figures  4-12  and  4-13  for  the  assumed  NPR  (7.0  and  16.0)  and  the  Mach  numbers  (1.9 
and  3.0).  For  supersonic  Mach  numbers,  the  pressure  distributions  show  the  characteristics 
of  the  pressure  rise  up  to  above  the  ambient  pressure  downstream  of  the  nozzle  throat  due 
to  the  separation  bubble  for  the  given  NPR. 

For  the  Mach  1.9  condition  and  NPR  7  0  (Figure  4-12),  the  exhaust  flow 
expansion  (lowest  to  highest  pressure  obtained)  downstream  of  the  nozzle  throat  for  the 
cowl  configuration  are  baseline,  -5  degree,  long,  and  -t-5  degree  cowl.  These  results  are 
dependent  on  the  internal  expansion  area  ratios  (the  nozzle  cowl  exit  divided  by  the 
nozzle  throat  area).  The  initial  compression  (lowest  to  highest  pressure)  exhibits  the  same 
pattern  as  the  initial  expansion.  The  flow  then  expands  and  recompresses.  The  flow 
reexpands  approximately  at  the  x=3.5  for  the  +5  degree  cowl.  jc=3.6  for  the  long  cowl, 
and  a=3.7  for  the  -5  degree  cowl  and  baseline  cowl,  respectively.  The  flow  recompresses 
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approximately  at  the  x=5.8  for  the  baseline  cowl,  ,t=6  2  for  the  long  cowl,  +5  degree 
cowl,  and  -5  degree  cowl,  respectively.  The  interaetion  between  the  boundary  layer  on  the 
top  surface  and  the  plume  shock  forms  the  other  separation  bubble  on  the  nozzle  surface. 
The  recompression  properties  depended  on  the  cowl  configuration  (nearly  in  the  reverse 
order  of  the  initial  compression)  which  decided  the  location  of  the  separation  bubble.  The 
influence  of  plume  shock  for  the  baseline  cowl  on  the  nozzle  surface,  i.e.,  the  location 
(x=5.8)  of  the  separation  bubble,  preceded  the  other  cases  u=6.2).  Therefore,  the  baseline 
cowl  recompresses  prior  to  the  cowl,  the  +5  degree  cowl,  and  the  -5  degree  cowl  on  the 
nozzle/surface 

Also,  for  the  Mach  3.0  and  NPR  15  0  condition  (Figure  4-13),  the  exhaust  flow 
expansion  t, lowest  to  highest  pressure)  downstream  of  the  throat  for  the  cowl  configuration 
are  baseline,  long,  and  +5  degree  cowl.  The  initial  compression  (lowest  to  highest 
pressure)  shows  the  same  sequence  as  the  initial  expansion.  After  compression,  the  flow 
reexpansion  properties  are  similar  and  result  in  nearly  the  same  (negative)  Cp.  For  Mach 
number  3.0,  the  cowl  shock  has  smaller  deflection  angle  than  Mach  number  1.9  so  that 
the  interaction  of  the  cowl  shock  with  boundary  layer  on  the  nozzle  surface  is  moved 
downstream  of  the  nozzle  surface.  Also  thickness  of  the  separation  bubble  is  thinner  than 
the  Mach  number  1.9  and  NPR  7  0.  As  s  result,  the  flow  appears  to  separation  without 
recompression. 

4.5  NPR  Effects 

NPR  effects  for  Mach  numbers  1.9  and  3.0,  for  the  baseline  cowl  configuration 


are  presented  in  Figures  4-14  and  4-15  Clearly,  the  supersonic  Mach  number 
characteristics  of  the  higher  pressure  of  the  initial  compression  for  the  increasing  nozzle 
pressure  ratios  arc  presented  in  Figures  4-11  and  4-12.  respectively. 

NPR  effects  for  Mach  number  1.9  (Figure  4-14)  have  the  properties  of  expansion, 
compression,  reexpansion,  and  recompression  for  all  cases.  The  apex  of  the  separation 
bubble  is  correspondingly  located  at  x=3.7  approximately  except  the  NPR=0.9,  i.e.,  flow 
off.  And.  the  thickness  of  the  separation  bubble  depends  on  the  NPR.  For  NPRs  7.0  and 
12.0,  the  exhaust  flow  expands  to  the  ambient  pressure.  As  NPR  increase,  the  location  of 
recompression  moves  downstream  of  the  nozzle  surface  from  ,r=3.9  (NPR=0.9)  to  x=l .2 
(NPR=12.))  due  to  the  length  of  the  separation  bubble. 

NPR  effects  for  Mach  number  3.0  (Figure  4-15)  exhibit  similar  trends  as  the  Mach 
number  1.9  near  the  initial  separation  bubble.  But.  as  mentioned  foi  the  nozzle/cowi 
effects,  the  cowl  shock  for  Mach  number  3.0  has  smaller  angle  than  Mach  number  1.9  so 
that  the  interaction  of  the  cowl  shock  with  the  boundary  layer  on  the  nozzle  surface  is 
moved  downstream  of  the  nozzle  surface.  Therefore,  the  exhaust  results  in  nearly  the 
same  (negative)  C„  and  appears  to  separate  without  recompression. 
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Figure  4-3  Convergence  History  (Explicit) 
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Figure  4-4  Grid  Resolution  Effects:  Long  Cowl,  M„=1.9,  NPR=7.0 
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Figure  4-5  Grid  Resolution  Effects:  Long  Cowl,  M„-3.0,  NPR=16 


Figure  4-6  Influence  of  Initial  Conditions:  Long  Cowl,  M„=1.9,  NPR=7.0 
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Figure  4-7  Mach  Number  Effects:  Baseline  Cowl,  NPR=l.O 
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Figure  4-8  Mach  Number  Effects:  Baseline  Cowl,  NPR=3.0 
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Figure  4-10  Mach  Number  Effects:  Baseline  Cowl,  NPR=7.0 
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Figure  4-11  Mach  Number  Effects:  Baseline  Cowl,  NPR=12.0 
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Figure  4-12  Nozzle  Cowl  Effects:  Baseline  Cowl,  M„=l.9,  NPR= 
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Figure  4-13  Nozzle  Cowl  Effects:  Baseline  Cowl,  M^=3.0,  NPR=16.0 
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Figure  4-14  NPB  Effects:  Baseline  Cowl.  M„=  1 .9 
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Figure  4-15  NPR  Effects:  Baseline  Cowl,  M„=3.0 
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V.  EXPERIMENTAL  AND  NUMERICAL  ANALYSIS 


For  the  nozzie/afterbodv  configuration  tested  through  experimental  and  numerical 
analysis,  two  .separation  bubbles  were  detected.  These  were  the  key  phenomena  to  be 
analyzed.  One  is  located  just  downstream  of  the  throat  and  is  due  to  a  sharp  corner  of  the 
nozzle'  throat,  while  the  other  is  located  on  the  top  surface,  where  an  interaction  of  ihe 
cowl  shocK  with  the  boundary  layer  of  the  top  surface  exwts.  Figure  5-1  is  a  diagram  of 
the  overexpanded  nozzle  exhaust  How  at  supersonic'  Mach  numbers.  According  to  the 
experimental  tests  t4).  the  oblique  shock  lor  Mach  number  1.9  and  3.0.  varies  from  45 
degrees  to  35  degrees,  respectively.  For  Mach  number  3  0.  the  cowl  shock  has  a  smaller 
deflection  angle  than  Mach  number  1  9  so  that  the  interaction  of  the  cowl  shock  with 
boundary  layer  on  the  nozzle  „urface  is  moved  downstream  of  the  nozzie  surface.  On  the 
assumption  of  the  computational  model,  the  exhaust  How  for  Mach  number  !  9  showed 
the  recompression  phenomena  at  the  second  separation  bubble.  But  the  How  for  Mach 
number  3.0  appeared  to  be  nearly  the  same  { negative)  C  and  separate  without 
recompression.  This  angle  changes  considerably  with  cowl  configuration  and  NPR . 

Unfortunately,  'here  were  pooi  agreement  between  the  experimental  data  and  the 
numerical  results  for  the  subsonic  cases.  These  results  were  caused  by  the  discrepancy 
between  the  models.  In  experimental  implementations.  Cochran  obtained  the  nozzle, 
surface  centerline  properties  for  the  performance  of  the  nozzle/afterbody  on  the 
assumption  of  the  two-dimensional  model,  although  he  did  the  three-dimensional  model. 
In  other  words,  his  results  included  the  three-dimensional  effects,  which  made  worse  the 
consistency  for  the  subsonic  cases  due  to  interaction  between  the  exhaust  nozzle  plume 
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and  the  external  flow.  And  lor  the  supersonic  eases,  especially  Mach  number  1.9,  due  to 
discrepancy  on  the  assumption  of  models,  the  three-dimensional  effects  can’t  be  neglected 
near  the  second  separation  bubble. 


This  angle  changes  considerably  with  cowi  configuration  and  NPR. 

5.1  Mach  Number  Effects 

5.1.1  Supersonic  Effects 

The  pressure  coefficient.  C?t  comparison  between  the  experimental  and  numerical 
method  were  presented  in  the  Figures  5-2  to  5-1 1  at  specific  NPR  for  the  baseline  cowl 
configuration.  Generally,  these  figures  show  a  favorable  agreement  between  theory  and 
experiment,  as  Mach  numbers  increased. 

For  the  NPR  1.0.  Figures  5-2  and  5-3  show  that  computational  versus  experimental 
results  are  a  little  more  consistent  for  Mach  number  3.0  than  for  Mach  number  1.9. 
Namely,,  for  the  Mach  number  1.9  (Figure  5-2)  according  to  the  numerical  result,  the 
pressure  distributions  appear  nearly  Hat  after  compression,  but  the  flow  recovers 
(recompresses  to  near  ambient  conditions;  at  the  last  pressure  tap  location  in  the 
experimental  results.  The  differences  between  the  experiment  and  numerical  resuits  are 
a  result  of  the  size  of  the  second  separation  bubble.  Foi  Mach  number  3.0,  the  numerical 
results  show  the  identical  trends  as  above,  but  the  experimental  results  have  the 
characteristics  of  slight  compression  and  reexpansion  around  the  second  separation  bubble 
U--6.0).  The  oblique  shock  angle  changes  considerably  with  Mach  number  so  that  Mach 


number  1.9  has  a  larger  reflection  angle  than  Mach  number  3.0.  This  reflection  angle 
determine  the  location  of  the  second  separation  area.  Comparing  Figures  5-2  and  5-3,  the 
location  of  the  recompression  for  Mach  number  1.9  is  ahead  of  that  of  Mach  number  3.0 
(.t=4.2  for  Mach  1.9  and  x=6.0  for  Mach  3.0). 

For  the  NPRs  3  0  and  5  0.  Figure  5-4  to  5-7  show  similar  trends  as  the  NPR  1.0. 
Figure  5-4  for  Mach  number  1  9  *  NPR  3  0  >.  on  the  assumption  if  the  laminar  flow,,  shows 
some  compression  characteristics  which  were  not  apparent  in  the  expenme  nal  results  near 
the  first  separation  bubble.  But,  the  experimental  results  show  a  higher  rate  of 
compression  around  the  second  separation  bubble  which  is  built  up,  contrary  fo 
expectation,  by  the  external  flow  (three-dimensional  effects)  as  well  as  the  interaction  of 
nozzle  surface  boundary  layer  and  cowl  shock  wave.  For  the  NPR  5  0,  the  difference  in 
\alues  of  the  pressure  coefficients  near  the  second  separation  bubble  as  well  as  near  the 
first  bubble  were  larger  than  the  NPR  3.0  cases.  Particularly,  at  x  values  beyond  the 
recompression,  the  difference  between  experiment  and  the  computational  solutions  grow 
larger.  The  Mach  number  3.0  results  show  a  better  agreement  than  those  at  Mach  1.9. 
Especially,  the  trend  of  reexpansion  and  recompression  after  the  apex  of  the  second 
separation  bubble  was  examined  at  the  NPR  5.0  (Figure  5-7)  for  the  experimental  results, 
because  the  external  flow  influenced  the  second  separation  bubble,  so  that  the  compressed 
up  to  x=5.9  and  reexpanded  slightly  alter  the  apex  of  the  second  separation  bubble.  This 
was  not  detected  in  the  computational  implementations  which  was  assumed  on  the  two- 
dimensional  model. 

For  .he  NPR  7  0  and  12.0,  Figures  5-8  to  5-li  show  that  the  trends  proceed 


72 


favorably  as  the  Mach  number  Increased,  although  the  numerical  investigation  was  based 
on  the  laminar  flow  assumption.  In  the  experimental  performance,  the  internal  flow 
initially  dominated  the  external  flow  downstream  of  the  nozzle  throat  for  the  same  NPR 
conditions,  as  Mach  number  increased.  Therefore,  there  was  better  agreement  for  Mach 
number  3.0  (Figures  5-0  and  5- ID  than  for  Mach  number  1.9  (Figures  5-8  and  5-10)  at 
the  second  separation  Kubble  as  well  as  the  initial  separation  bubble.  Additionally,  these 
results  mean  that  the  difference  between  the  turbulence  and  the  laminar  flow  make  worse 
the  results  between  experiment  and  theory  at  the  initial  separation  bubble  as  the  Mach 
number  decreases.  Also,  the  difference  between  the  two-dimensional  model  and  the  three- 
dimensional  model  increases  near  the  second  separation  bubble,  since  the  internal  flow 
is  effected  by  the  external  flow  so  that  the  exhaust  How  -  the  results  were  extracted  along 
the  centerline  of  tne  nozzle  surface  at  the  experimental  performance  -  continuously 
compressed  at  this  region.  Also,  as  the  NPR  increased,  the  agreement  between  theory  and 
experiment  priced  unfavorable  locally  for  the  same  Mach  number  conditions  at  the  apex 
of  the  initial  separation  bubble,  since  the  flow  with  the  higher  NPR  formed  the  larger 
separation  bubble  at  this  region. 

To  investigate  tne  difference  between  the  experiment  and  the  numerical  solutions 
in  detail,  the  three-dimensional  effects  were  considered  for  specific  cases.  Figures  5-12 
and  5- 13  show  the  comparison  between  the  experimental  results  and  the  numerical  results. 
The  experimental  resuits  were  extracted  from  the  centerline  and  off-centerline  (positive 
1  inch  and  negative  1  inch  laterally)  tap  location  on  the  nozzle  surface.For  Mach  number 
1.9  and  NPR  3.0,  Figure  5-12.  the  internal  How  initially  dominates  the  external  flow 


downstream  of  the  nozzle  throat  so  that  off-centerline  pressure  above  the  centerline  one. 
But  after  approximate!)  *=4.8,  the  centerline  flow  maintained  the  higher  pressure  than  the 
off-center  line  and  compressed  by  the  influence  of  the  external  flow,  But,  the  off- 
centerline  flow  slightly  compressed  and  expanded  after  *=4.7.  and  the  flow  continuously 
compressed  again  due  to  the  adverse  flow  on  the  nozzle  surface  at  *=6.2.  Mach  number 
3.0  and  NPR  3.0  case.  Figure  5-13.  shows  that  there  are  smaller  variations  between  the 
experiment  and  the  numerical  solutions,  because  the  internal  flow  dominates  the  external 
flow  as  Mach  number  increases.  As  a  result,  as  Mach  number  increased,  the  conversion 
point  where  the  signal  of  the  difference  value  between  the  centerline  and  the  off-centerline 
is  changed  moves  afterward  (from  ,c=4  7  to  *=5.5),  which  explains  why  the  internal  flow 
dominates  the  external  flow  and  there  is  a  thin  boundary  layer  on  the  nozzle  surface. 

5.1.2  Subsonic  Effects 

The  subsonic  cases  showed  poor  agreement  between  the  experiment  and  the 
computational  solutions,  because  there  is  severe  interaction  between  the  nozzle  exhaust 
plume  and  the  external  flow.  To  investigate  the  difference  between  the  experiment  and 
the  computational  tesults.  the  three  dimensional  effects  were  considered  for  specific  case 
like  the  supersonic  ones.  For  Mach  number  0.6  and  NPR  3.0,  Figure  5-14  shows  the 
comparison  between  the  experiment  and  the  computational  results.  Clearly,  the 
experimental  results  showed  the  complicated  flow  phenomena  on  the  nozzle  surface 
comparing  to  the  results  for  the  supersonic  cases  which  showed  the  symmetric  properties 
for  the  off-centerline  values.  Just  downstream  of  the  nozzle  throat,  the  internal  flow  is 
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influenced  by  the  external  so  that  centerline  pressure  is  higher  than  the  off-centerline 
pressure.  And  then  the  centerline  pressure  ts  lower  than  the  off-centerline  one,  which 
means  that  initially  the  flow  turns  towards  the  centerline  and  then  turn  outward  after  the 
impingement  at  the  centerline.  The  centerline  flow  promptly  recovered  the  higher  pressure 
than  the  off-centerline  and  maintained  nearly  above  the  ambient  pressure  uniformly,  after 
ail.  There  is  a  large  difference  between  the  results  at  the  initial  separation  bubble,  which 
means  that  the  assumption  of  turbulent  versus  laminar  flow  is  more  significant  in  the 
subsonic  regime.  To  reinforce  the  subsonic  problems,  some  supplementary  methods  were 
used,  w  hich  are  the  assumption  of  the  lower  (external)  fictitious  boundary  conditions. 
These  methods  still  could  not  solve  the  subsonic  problems. 

5.2  Nozzle  Cowl  Effects 

The  pressure  coefficient,  Cp ,  comparison  between  the  experimental  and  numerical 
methods  are  presented  in  the  Figures  5-15  to  5-21  at  specific  NPRs  for  the  different  cowl 
configuration. 

The  case  for  the  long  cowl  at  Mach  number  1.9  and  NPR  7.0,  Figure  5-16,  shows 
that  the  differences  in  pressure  coefficient  between  theory  and  experiment  around  the  apex 
of  the  compression  are  larger  than  the  baseline  cowl.  Figure  5-15.  This  is  because  the 
expansion  waves  caused  by  the  sharp  throat  impinge  on  the  upper  surface  of  the  cowl  and 
the  resulting  exhaust  flow  is  much  more  complicated.  Therefore,  the  numerical  results 
show  larger  separation  bubble  than  the  experimental  test.  Since  a  turbulent  boundary  layer 
allows  higher  adverse  flow  gradient,  there  is  smaller  separation  bubble  on  the  downstream 
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side  of  the  nozzle  surface.  Also,  the  identical  results  were  obtained  for  the  +5  degree  cowl 
in  the  Figure  5-17.  For  the  -5  degree  cowl.  Figure  5-18,  the  exhaust  flow  is  expanded  at 
the  hinge  point  of  the  cowl  upper  surface,  so  that  there  will  be  smaller  separation  bubble 
than  the  long  and  +5  degree  cowl.  Therefore,  the  results  for  -5  degree  cowl  case  show  a 
better  agreement  than  the  long  and  +5  degree  cowl. 

The  case  for  Mach  number  3.0  and  NPR  16.0  shows  similar  trends  as  above  until 
the  flow  encounters  compression,  and  it  again  shows  good  agreement  after  reexpansion. 
In  contrast,  for  Mach  number  1 .9  and  NPR  7.0  case,  the  numerical  results  appear  Hat  after 
compression,  while  the  experimental  results  show  the  recompression.  As  the  Mach  number 
increased,  the  internal  flow  dominated  the  external  How  so  that  the  internal  flow  was  not 
influenced  by  the  external  flow  and  kept  the  uniform  flow. 

5.3  NPR  Effects 

The  pressure  coefficient.  Cp,  comparison  between  the  experimental  and  numerical 
method  were  already  presented  in  the  Figures  5-3  to  5-15,  Figures  5-19,  5-22,  5-23,  and 
5-24  for  various  NPR  for  the  baseline  cowl  configuration.  Generally,  it  showed  a  good 
agreement  between  the  numerical  solutions  and  experiment  as  Mach  number  and  NPR 
increased. 

For  Mach  number  1 .9,  the  Figures  5-2,  5-4,  5-6,  and  5- 10  show  that  the  trends  are 
less  consistent  than  for  Mach  3.0,  Figures  5-7,  5-9,  5-22,  5-11,  5-19,  5-23,  and  5-24. 
However,  there  are  some  discrepancies  between  theory  and  experiment  around  the  peak 
of  the  compression  and  the  second  separation  bubble  exist,  because  the  boundary  layer. 


most  likely  turbulent  in  the  experiment,  allow  a  higher  adverse  flow  gradient.  In  addition, 
it  should  be  neeessary  to  consider  the  other  factors  that  effect  the  phenomenon,  i.e„  the 
disagreement  at  x  values  beyond  the  recompression.  The  detailed  explanation  to 
synthesize  the  whole  results  will  be  covered  below.  Figures  5-7.  5-9,  5-22,  5-1 1,  5-19, 
5-23,  and  5-24,  for  the  Mach  number  3.0,  show  that  the  agreement  between  numerical 
soiuuuns  and  experiment  proceed  quite  favorably,  as  the  NPR  increases,  although  the 
numerical  investigation  was  performed  on  the  basis  of  a  laminar  How  assumption. 

In  the  experimental  testing,  Cochran  designed  the  three-dimensional  model  and 
conducted  the  test.  However,  he  assumed  a  two-dimensional  hypersonic  nozzle /afterbody 
for  the  performance  parameters  of  the  test,  and  reduced  data  according  to  the  parametric 
cases.  Even  though  he  fully  performed  the  three-dimensional  model  ,  he  calculated  the 
centerline  of  integrated  external  ramp  pressure  coefficient.  Cp,  on  the  assumption  of  two- 
dimensional  model,  since  the  higher  density  of  pressure  taps  characterized  the  nozzle 
exhaust  flow  and  the  afterbody  ramp  pressure  more  accurately.  Additionally,  he  obtained 
the  off-centerline  pressure  to  determine  the  lateral  variation  of  the  flow,  which  will  be 
valuable  data  for  the  further  research.  Therefore,  the  three-dimensional  effects  were 
performed  to  validate  the  consistency  between  the  experimental  data  and  the  numerical 
solutions.  Generally,  as  Mach  number  increased,  the  difference  between  the  experiment 
and  the  numerical  results  decreased  favorably.  For  the  subsonic  case,  there  was  severe 
interaction  between  the  nozzle  plume  flow,  the  external  How,  and  the  afterbody  flow.  This 
interaction  made  comparison  of  the  CFD  data  and  experimental  data  very  difficult,  so  that 
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the  comparison  for  just  on  case  was  conducted.  But,  for  the  supersonic  cases,  the  results 
show  good  agreement  near  the  initial  separation  bubble  as  Mach  number  increased, 
because  the  exhaust  flow  overexpansion  dominates  and  is  uninfluenced  by  other  effects. 
Near  the  second  separation  bubble,  however,  the  results  are  not  as  consistent,  because  the 
recompression  characteristics  are  intiuenced  by  the  interaction  among  the  external  flow 
all  around,  the  afterbody  plume,  as  well  as  the  separation  bubble. 

For  the  Mach  number  effects,  as  Mach  number  increased,  there  is  better  agreement 
between  the  experimental  and  the  numerical  results  at  the  initial  separation  bubble  than 
at  the  second  separation  bubble,  since  the  internal  properties  are  initially  uninfluenced  by 
the  external  flow.  Additionally,  the  assumption  of  the  laminar  flow  slightly  caused  the 
differences  between  the  data.  Also,  for  Mach  number  1.9  cases,  the  discrepancy  between 
the  two-dimensional  computational  model  and  the  three-dimensional  experimental  model 
caused  the  difference.  For  the  nozzle/cowl  configuration,  the  -5  degree  cowl  showed  the 
best  agreement  among  the  cowl  configurations,  because  the  downward  cowl  (5  degrees) 
produced  a  more  dominant  exhaust  flow  downstream  of  the  nozzle  throat  and  more 
effectively  isolated  the  nozzle  from  its  surroundings.As  for  the  Mach  number  effects,  the 
laminar  flow  assumption  in  the  CFD  model  had  its  largest  influence  at  low  Mach 
numbers,  with  decreasing  influence  as  Mach  number  increased.  But  for  the  NPR  effects, 
the  flow  with  higher  NPR  formed  the  thicker  initial  separation  bubble. 

If  it  is  assumed  that  the  research  was  performed  on  the  same  two-dimensional 
model,  we  may  analyze  and  compare  the  results  more  accurately.  As  stated  in  reference 
(18),  the  laminar  boundary  layer  on  the  top  surface  separates  because  of  the  interaction 
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of  the  cowl  shock  with  the  boundary  layer  on  the  top  surface.  But  separation  completely 
disappeared  for  the  turbulent  flow  under  the  experimental  test.  The  turbulent  boundary 
layer  is  able  to  accept  a  higher  adverse  flow  gradient  without  separating  for  the  two- 
dimensional  assumption  model.  Therefore,  the  numerical  investigation  for  the  three- 
dimensional  model  can  solve  the  obscure  problems  encountered  near  the  secondary 
separation  bubble. 
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Figure  5-2  Comparison  of  an  Experimental  and  Numerical  Pressure  Coefficients 
(Mach  Effects):  Baseline  Cowl,  M„=1.9,  NPR=l.O 
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Figure  5-3  Comparison  of  an  Experimental  and  Numerical  Pressure  Coefficients 


(Mach  Effects):  Baseline  Cowl,  M„=3.0,  NPR=1.0 
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Figure  5-4  Comparison  of  an  Experimental  and  Numerical  Pressure  Coefficients 
(Mach  Effects):  Baseline  Cowl,  M„=  1 .9,  NPR=3.0 
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Comparison  of  an  Experimental  and  Numerical  Pressure  Coefficients 
(Mach  Effects):  Baseline  Cowl,  M„=3.0,  NPR=3.0 
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Figure  5-6  Comparison  of  an  Experimental  and  Numerical  Pressure  Coefficients 


(Mach  Effects):  Baseline  Cowl,  M„=1.9,  NPR=5.0 
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Figure  5-7  Comparison  of  an  Experimental  and  Numerical  Pressure  Coefficients 
(Mach  Effects):  Baseline  Cowl,  M„=3.0,  NPR=5.0 
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Figure  5*8  Comparison  of  an  Experimental  and  Numerical  Pressure  Coefficients 
(Mach  Effects):  Baseline  Cowl,  M^=1.9,  NPR=7.0 
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Figure  5-10  Comparison  of  an  Experimental  and  Numerical  Pressure  Coefficients 
(Mach  Effects):  Baseline  Cowl,  M„=1.9,  NPR=12.0 
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Figure  5-11  Comparison  of  an  Experimental  and  Numerical  Pressure  Coefficients 
(Mach  Effects):  Baseline  Cowl,  M„=3.0,  NPR=12.0 
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Figure  5-12  Three-Dimensional  Effects  for  the  Supersonic  Gases:  Baseline  Gowl, 
Mj=1.9,  NPR=3.0 
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Figure  5-13  Three-Dimensional  Effects  for  the  Supersonic  Gases:  Baseline  Cowl, 
M„=3.0,  NPR=3.0 
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Figure  5-14  Three-Dimensional  Effects  for  the  Subsonic  Cases:  Baseline  Cowl, 
M„=0.6,  NPR=3.0 
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Figure  5-15 


Comparison  of  an  Experimental  and  Numerical  Pressure  Coefficients 
(Nozzle  Cowl  Effects):  Baseline  Cowl,  M„=1.9,  NPR=7.0 
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Figure  5-16  Comparison  of  an  Experimental  and  Numerical  Pressure  Coefficients 
(Nozzle  Cowl  Effects):  Long  Cowl,  M„=1.9,  NPR=7.0 
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Figure  5-17  Comparison  of  an  Experimental  and  Numerical  Pressure  Coefficients 


(Nozzle  Cowl  Effects):  +5  Degree  Cowl,  M„=1.9,  NPR=7.0 
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Figure  5-18  Comparison  of  an  Experimental  and  Numerical  Pressure  Coefficients 
(Nozzle  Cowl  Effects):  -5  Degree  Cowl,  M„=1.9,  NPR=7.0 
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NOZZLE/COWL  PRESSURE  COEFFICIENT 

DASU1JNU  COWIiiTI,N184);  M.-3.0.  Nl’R»16.0;  RU.-2  03757  B06/IT 


Figure  5-19 


(Nozzle  Cowl  Effects):  Baseline  Cowl,  Mm=3.0,  NPR=16.0 


97 


Figure  5-20  Comparison  of  an  Experimental  and  Numerical  Pressure  Coefficients 
(Nozzle  Cowl  Effects):  Long  Cowl,  M,=3.0,  NPR=16.0 
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Mom-I.  STATION  IINCIIJ 


Figure  5-21  Comparison  of  an  Experimental  and  Numerical  Pressure  Coefficients 
(Nozzle  Cowl  Effects):  +5  Degree  Cowl,  M„=3.0,  NPR=16.0 
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MODEL  STATION  UNCH| 


Figure  5-22  Comparison  of  an  Experimental  and  Numerical  Pressure  Coefficients  (NPR 
Effects):  Baseline  Cowl,  M„=3.0,  NPR-9.0 
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(Nozzle  Cowl  Effects):  Baseline  Cowl,  M,=3.0,  NPR=16.0 
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Figure  5-24  Comparison  of  an  Experimental  and  Numerical  Pressure  Coefficients 
(Nozzle  Cowl  Effects):  Baseline  Cowl,  Mw=3.0,  NPR=1620.0 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 

6.1  Conclusions 

Off-design  flow  over  a  quasi  two-dimensional  hypersonic  nozzle  /afterbody  model 
was  computationally  solved  assuming  laminar,  planar  flow  to  validate  the  consistency  of 
the  experimental  data  by  Cochran  (4).  Two  flux-split  algorithms,  flux- vector  splitting  of 
Van  Leer  and  flux-difference  splitting  of  Roe,  were  explicitly  investigated  for  this 
research.  The  code  developed  by  WL/FIMC  performed  quite  well,  its  finite  volume 
formulation  enhancing  the  practical  stability  with  skewed  cell  geometries. 

Generation  of  the  two-dimensional  grid  was  a  major  effort,  made  simpler  by  using 
the  GRIDGEN  package.  The  elliptic  solver  was  used  to  maintain  smooth  grid  contours 
for  the  finite  volume  scheme.  In  order  to  provide  accurate  resolution  of  the  shock  region 
and  the  boundary  layer,  the  grid  points  were  geometrically  clustered. 

The  two-dimensional  explicit  flux  splitting  algorithms  performed  very  well.  The 
finite  volume  implementation,  which  the  code  used  for  an  extremely  robust  algorithm, 
handled  cells  with  very  small  volumes  as  well  as  regions  with  highly  skewed  cells.  The 
code  generally  overpredicted  the  experimental  data,  with  deviations  attributable  to 
numerical  error  in  either  the  afterbody  model  or  in  the  numerical  solution  procedure.  In 
particular,  the  laminar  assumption  used  for  the  computational  model  did  not  accurately 
capture  the  initial  separation  bubble  which  was  most  likely  turbulent  in  the  experimental 
tests.  Also,  the  two-dimensional  model  used  in  the  computational  solutions  showed 
inconsistency  at  the  secondary  separation  bubble  due  to  the  interaction  between  the  nozzle 
plume  and  the  external  flow,  i.e.,  the  three-dimensional  effects. 


As  expected,  the  nozzle/ afterbody  exhaust  flow  was  overexpanded  and  a  separation 
bubble  occurred  just  downstream  of  the  throat  due  to  the  sudden  increase  of  the  area  in 
the  divergent  ramp  and  the  sharp  comer  at  the  throat.  In  addition,  the  second  separation 
bubble  was  caused  by  the  interaction  of  the  cowl  shock  and  the  boundary  layer  of  the  top 
surface. 

The  flow  recompression  and  reexpansion  properties  over  all  of  the  nozzle 
/afterbody  depended  on  the  parametric  cases  and  the  grid  refinement.  Generally,  as  the 
Mach  number  and  NPR  (Nozzle  Pressure  Ratio)  increased,  there  was  increasing  agreement 
between  the  experimental  and  the  computational  results,  even  though  the  computational 
investigation  was  performed  oh  the  basis  of  the  laminar  flow  assumption.  For  the  +5 
degrees  and  long  cowl  configurations,  the  agreement  worsened,  because  of  the 
complicated  flow  structure  and  separation  region  near  the  sharp  throat  and  the  upper 
surface  of  the  cowl.  A  refined  grid  showed  that  it  is  possible  to  analyze  the  phenomena 
of  the  flow  in  detail  and  capture  the  scale  of  the  separation  bubble  more  accurately  near 
the  downstream  of  the  throat.  The  discrepancy  of  methodology  and  assumption  between 
the  experiment  and  the  theory  explains  the  observed  variation  in  the  results. 

6.2  Recommendations 

This  thesis  will  provide  a  basis  for  the  solution  of  a  specified  nozzle/afterbody 
problems  in  further  research.  The  following  recommendations  are  the  result  of  careful 
consideration  of  the  research  and  conclusions  reached  in  the  course  of  this  thesis. 

In  Cochran's  experimental  testing,  a  three-dimensional  hypersonic  nozzle  /afterbody 
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model  was  investigated,  but  the  pressure  coefficients  for  the  performance  parameter  were 
calculated  oh  the  assumption  of  a  two-dimensional  model.  Also,  this  numerical  research 
was  computationally  conducted  assuming  a  two-dimensional  model.  Some  discrepancies 
between  Cochran's  three-dimensional  and  the  current  two-dimensional  models  brought 
about  some  inconsistency  and  require  further  research.  Therefore,  it  is  desirable  to 
investigate  the  three-dimensional  computational  model  and  validate  the  experimental 
results.  Furthermore,  future  plans  should  include  turbulence  modeling  and  validation  of 
the  revised  code  with  nozzle/afterbody  experimental  data. 
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APPENDIX  A:  CN1DAT  CONTROL  FILE 
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IL,  JL,  IMETRC 
102  72  0 

NEND 

1000 

ICFL,  CFLEXP,  CFLMAX,  CFL 
5  1  0.9  0.01 

IREST,  CFCRHO,  CFCEI,  CFLPEN,  CEXPPEN,  INOFRZ 
1  5.00  5.00  1  1  1 

IMPLTJLGTSTJSWVLJLMTR,OMEGAX)ELTEPX>ELTILJENTHJNS^SWPS 
0  14  2  1  l.E-10  0.5  2  1  2 

IADBWL,  ICASE,  TWALL,  ALPHA,  RM,  REL,  RL,  TINF,  IGRID 
0  2  400.6  0.  1.9026  2.00984E6  1.0  105.192  2 

IREAD,  IP3DOP,  IDGBUG,  MODPR,  NRST,  EFMRTI,  IFMRTO,  IINT1 
1  2  0  5  0  1  1  1 

INS1,  INS2,  JNS1,  JNS2 
1  52  37  46 

TINWALL,  AINLPHA,  RINM,  RINEL,  RINL,  TININF  IGRID 
400.6  0.  0.29892  1.25155E+6  1.0  481.56  2 
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APPENDIX  B:  NOZZLE/AFTERBODY  DATA 


no 


able  B-l  Mach  Number  Effects 


V". 


*  TPN  :  Test  Point  Number,  PN  :  Part  Number 

*  Subcript  oo  =  External  Freest!  cam  Condition,  n  =  Internal(Noz/.le)  Freeslream  Conditon 


Nozzle  Cowl  Effects 


Table  B-3  NPk  Effects 


2. 09480  l  8.99580  2.00542  53.717  200.78  0.28331  0.47941  457.030  480.35 


APPENDIX  C:  PRESSURE  CONTOUR 


Figure  C-l  Nozzle/Afterbodv  Pressure  Contour:  Baseline  Cowl,  Mw=1.9,  NPR=1.0 
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Noz/le/Al'terbody  Prcs.su rc  Contour:  Baseline  Cowl,  M„=3.0,  NPR=1.0 
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Figure  C-4  Nozzle/Afterbody  Pressure  Contour.  Baseline  Cowl,  M„=3.0,  NPR=3.0 
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Figure  C*5 
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Nozzle/Afterbody  Pressure  Coniour:  Baseline  Cowl,  M„=1.9,  NPR=5.0 
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Figure  C-6  Nozzle/ Afterbody  Pressure  Contour.  Baseline  Cowl,  M^=3.0,  NPR=5  0 
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Figure  C-7  Nozzle/Al'terboUv  Pressure  Contour  Baseline  Cowl,  M„=1.9,  NPR=7.0 
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Figure  C-8  Nozzle/Afterbody  Pressure  Contour:  Baseline  Cowl,  M„=3.0,  NPR=7.0 
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Figure  C-9  Nozzle/Afterbody  Pressure  Contour:  Baseline  Cowl,  M„=1.9,  NPR=12.0 
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Figure  C-10  Nozzle/Afterbody  Pressure  Contour.  Baseline  Cowl,  M„=3.0,  NPR=12.0 


X  [INCH) 


Figure  C-14  Nozzle/Afteibody  Prcj.su re  Contour.  Baseline  Cowl,  Mw=3.0,  NPR=16.0 
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Figure  C-15  Nozzle/Afterbodv  Pressure  Contour:  Long  Cowl,  M,=3.0,  NPR=16.0 
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APPENDIX  D:  MACH  NUMBER  CONTOUR 
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Figure  D-2  Nozzle/ Afterbody  Mach  Number  Contour:  Baseline  Cowl,  M^=3,0,  NPR-1.0 
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Figure  D-4  Nozzle/Afterbody  Mach  Number  Contour:  Baseline  Cowl,  M^=3.0,  NPR=3.0 
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Figure  D-5  Noz/ie/Alterbody  Mach  Number  Contour:  Baseline  Cowl,  M„=1.9,  NPR=5.0 
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NOZZEK/COWE .MACII  NUMBER  CONTOUR 


Figure  D-8  Nozzle/ Afterbody  Mach  Number  Contour:  Baseline  Cowi, 
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Figure  D-ll  Nozzle/Afterbody  Mach  Number  Contour:  Long  Cowl,  M„=1.9,  NPR=7.0 
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Figure  D-12  Nozzle/Afterbody  Mach  Number  Contour:  +5  Degree  Cowl,  M„=1.9,  NPR=7.0 
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Figure  D-15  Nozzle/Afterbody  Mach  Number  Contour:  Long  Cowl,  Mj=3.0,  NPR=16.0 
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This  computational  study  of  a  hypersonic  sub-scale  nozzle/cowl  config¬ 
uration  compares  numerical  solutions  with  experimental  data  for  nozzle 
performance  at  off  design  conditions  (subsonic  through  supersonic  Mach 
numbers.  A  combination  of  Van  Leer  and  Roe  flux-splitting  algorithms 
was  used  to  solve  the  flow  field.  The  explicit  formulation  assumed 
laminar  planar  flow  and  the  perfect  gas  equation  of  state.  A  valida¬ 
tion  for  the  numerical  algorithm  was  accomplished  by  comparing  with 
experimental  data  for  a  variety  of  cases.  Grid  refinement,  geometry, 
and  operating  conditions  all  influenced  the  comparison  of  the  experi¬ 
mental  and  computational  data.  Generally,  as  the  Mach  number  and 
nozzle  pressure  ratio  increased,  the  consistency  between  the  experimen¬ 
tal  results  showed  better  agreement,  although  the  computational  inves¬ 
tigation  was  performed  on  the  basis  of  planar,  laminar  flow.  For  the 
different  cowl  configurations,  the  results  for  baseline  and  -5  degree 
cowl  show  a  better  agreement  than  the  long  and  +5  degree  cowl. 
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